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ABSTRACT 

In this work, we present a comprehensive X-ray picture of the interaction between a super star 
cluster and the ISM. In order to do that, we compare and combine the X-ray emission from the 
superwind driven by the cluster with the emission from the wind-blown bubble. Detailed ana- 
lytical models for the hydrodynamics and X-ray luminosity of fast polytropic superwinds are 
presented. The superwind X-ray luminosity models are an extension of the results obtained 
in Paper I. Here, the superwind polytropic character allows to parameterize a wide variety 
of effects, for instance, radiative cooling. Additionally, X-ray properties that are valid for all 
bubble models taking thermal evaporation into account are derived. The final X-ray picture is 
obtained by calculating analytically the expected surface brightness and weighted temperature 
of each component. All of our X-ray models have an explicit dependence on metallicity and 
admit general emissivities as functions of the hydrodynamical variables. We consider a real- 
istic X-ray emissivity that separates the contributions from hydrogen and metals. The paper 
ends with a comparison of the models with observational data. 
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1 INTRODUCTION 

Young super star clusters (SSCs) are compact stellar aggregates resulting from intense, concentrated, violent star formation episodes. With 
ages between ~ 10 6 -10 7 yr and masses in the range of ~ 10 4 -10 7 Mq, they represent the current analogues of globular clusters and 
have been identified as the building blocks of stellar formation in many galaxies (O'connell et al. 1995; Whitmore & Schweizer 1995; Ho 
1997; Whitmore et al. 1999; Melo et al. 2005 and references therein). Within its effective radius of up to a few parsecs (< 10-pc, Meurer 
et al. 1995; Melo et al. 2005), a SSC contains many thousands or tens of thousands of early-type stars (Leitherer & Heckman 1995) that 
dominate the injection of mass and energy (~ 10 38 -10 41 erg s _1 ) to the ISM through the overall contribution of their individual stellar 
winds and supemovae explosions. The interactions that take place among these inner flows lead to the thermalization of their kinetic energy 
through strong shock waves that generate a high central over-pressure. As a result, a powerful superwind is driven out of the cluster at large 
speed (> 1000 km s _1 ). Outside of the cluster, the superwind expands almost freely until the presence of the ISM becomes important. 
The subsequent interaction of this supersonic and metal-rich outflow with the ambient interstellar gas generates a structure with a more 
complex dynamics, namely, a superbubble (Castor, McCray & Weaver 1975; Weaver et al. 1977, W77 hereafter; Koo & McKee 1992a,b; 
Bisnovatyi-Kogan & Silich 1995 and references therein). The configuration of a superbubble is traced by two shock fronts that are inherent 
to its hydrodynamical evolution: a leading (outer) shock that sweeps up, compresses and confines the ambient ISM to an external shell; and 
a secondary (inner) shock that propagates backwards (in the leading shock reference frame) and thermalises the kinetic energy of the free 
superwind. The resulting shocked gases are separated by a contact discontinuity and the thermal pressure of the hot, shocked superwind 
drives the outer shell. Soon, the mass of the shocked ambient gas becomes progressively larger than that of the shocked superwind. As this 
tendency continues, the outer shock becomes radiative, and the shocked ISM collapses into a thin, dense shell. Being more tenuous, the 
bubble hot interior has a larger cooling time scale and thus it remains quasi-adiabatic. This stage is called the snowplow phase. Here, in the 
same fashion as in Paper I, the resulting structure is split into the four zones shown in Fig.fT] region A corresponds to the star cluster volume 
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Figure 1. Structure of a superbubble. Each zone is defined according to its contents: region A contains the thermalised SNe and stellar wind ejecta, region B 
contains the free-expanding wind, region C represents the hot superbubble where the wind thermalised at the inner shock is confined, and region D is the shell 
that contains the swept-up ambient gas. An unperturbed ISM surrounds the structure. 

itself, region B is the zone where the wind freely expands, region C is the superbubble that contains the shocked superwind, and region D is 
the outer shell where the leading shock accumulates ambient gas. 

Due to the high temperatures produced by the multiple shock waves and the large amount of mass deposited by the sources, both 
superwinds and superbubbles are potential diffuse X-ray emitters (Chevalier 1992; Chu et al. 1995, C95 henceforth; Canto, Raga & Rodriguez 
2000; Stevens & Hartwell 2003). This represents an opportunity to study the X-ray emission associated with massive stellar clusters to better 
understand and constrain the ongoing hydrodynamics and its energetics. With the aid of Chandra, it has been found that several clusters in 
the Local Group emit an important amount of diffuse X-rays: near the centre of our Galaxy, the X-ray emissions of the Arches and Quintuplet 
clusters have been detected and studied by Yusef-Zadeh et al. (2002) and Law & Yusef-Zadeh (2004); that of NGC 3603 -one of the most 
luminous clusters in the Milky Way- has been discussed by Moffat et al. (2002); and Stevens & Hartwell (2003) have analysed R136 and 
NGC 346 in the Magellanic Clouds (LMC and SMC, respectively). On the other hand, Chu & Mac Low (1990) examined the Einstein 
archives for LMC and found seven superbubbles that were related to OB associations. Repeating later the same analysis using ROSAT data, 
they found four additional but dimmer superbubbles (C95). More recently, Townsley et al. (2003) have reported diffuse emission from a 
superbubble associated to M 17, also known as the Omega Nebula; similarly, Smith & Wang (2004) have used XMM-Newton data to study 
the emission from a 100-pc superbubble related to the dense OB association LH90 in 30 Dorados. In distant galaxies the X-ray emission 
presents two components: a system of point-like non-resolved sources, and an extended diffuse component possibly associated to a collection 
of superbubbles generated by individual SSCs (Summers et al. 2004; Smith, Struck & Nowak 2005). 

This theoretical work has a twofold objective. The first one is to present detailed hydrodynamic and X-ray luminosity models for fast 
polytropic superwinds, models that can be easily modified to account for general emissivity functions (Section[2]l. This property allowed us 
to explicitly separate the X-ray contributions from hydrogen and metals and to obtain the plasma weighted temperature. These results could 
be of particular use when comparing with synthetic or observed spectra, since we obtained them theoretically taking into consideration bands 
suitable to the instruments on board Chandra and XMM-Newton. Additionally, we also present analytical formulas for the corresponding 
X-ray surface brightness. The second objective is to undertake the task of obtaining thoroughly all of the above for the standard analytical 
bubble models (Section|3j- Three properties that are universal for bubbles accounting for shell evaporation are derived: an X-ray luminosity 
scaling factor that is independent of the particular dynamical evolution, the shape of the related surface brightness profile, and the projected 
X-ray temperature. We examine also the case of when the evaporation of material from the outer shell is prohibited (Section [4]l. Later, we 
assemble our models and present the complete X-ray panorama (Sections[5}. Finally, we make a comparison with observational data (Section 

All of our X-ray models are analytical. Moreover, they are independent of the actual shape of the emissivity function because they involve 
hypergeometric functions. In the more general case, they can be evaluated using the code of Colavecchia & Gasaneo (2004); however, that 
would not be necessary here because we have constructed closed-form expansions in terms of elementary algebraic expressions. From the 
theoretical point of view, the freedom of selecting the emissivity function and the separation of the contributions from hydrogen and metals 
are advantages, since it is usual in the literature to assume a constant X-ray emissivity (C95; Garcfa-Segura & Mac Low 1995) and that the 
X-ray luminosity scales linearly with metallicity (oc Z). The first advantage enables a more realistic modeling, and the second allows to 
discriminate the importance of each contribution at a given temperature and metal abundance. Our formalism also permits a straightforward 
calculation of other integral properties such as the luminosity of the plasma in other bands and the mass of the emitting gas (see Appendix 
IA2b . On the hydrodynamical side, the polytropic character of the winds makes it possible to implicitly parameterize the effect of radiative 
cooling in a simple manner. 

Throughout the Paper, astrophysical units are used when convenient: ta is the time in Myr, Es$ is the mechanical luminosity in units of 
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10 38 erg s^ 1 , i? sc ,pc is the star cluster radius in parsecs, V$ is the terminal speed in units of 1000 km s -1 , Tk c v is the plasma temperature at 
the star cluster centre in units of keV and Ax, -23 is the X-ray emissivity in units of 10 -23 erg s _1 cm 3 . A comprehensive list of symbols is 
given in Table ID ll The reader not interested in the detailed theoretical derivations may wish to proceed directly to Sections[5]and[6] 



2 POLYTROPIC SUPERWINDS 

Stationary superwind models have proven valuable in reproducing the average behaviour of the star-cluster outflows obtained from full 3D 
hydrodynamical simulations. Although the continuous wind-source distributions assumed in the analytical models do not produce the detailed 
features that might result from the local interaction of stellar winds and supernovae explosions, they provide a good approximation for the 
radially averaged hydrodynamical profiles even for clusters with only ~100 stars and become an increasingly accurate approximation as the 
cluster mass and stellar density also increase. Furthermore, not to far away from the cluster the agreement is better than in the inner zone 
because the inhomogeneities introduced by the discrete stellar population are quickly smoothed out (Canto et al. 2000). 

After a short initial transitory period, both the averaged hydrodynamical profiles and the X-ray luminosity reach almost steady values. 
The transitory period lasts until the cluster core is filled with shocked gas, which depends on the cluster radius and mass and energy injection 
rates. For instance, in their simulations, Rockefeller et al. (2005) have found this time to be ~ 2 x 10 4 yr for the Arches cluster (R BC , P c ~ 0.2 
pc), whereas for the ~ 5 times larger Quintuplet cluster they found it to be larger than ~ 10 5 yr. They also found small temporal variations 
of the X-ray luminosity: ~ 1 percent for the Arches cluster and around 4-7 per cent for the Quintuplet cluster over a period of ~ 50 yr. 

All of the above implies that stationary models can be used to sequentially approximate the hydrodynamical evolution of a cluster as 
well as the corresponding diffuse X-ray emission, given that the self-adjustment times between transitions are short and not too frequent in 
comparison with the cluster age and the intrinsic temporal variations of the quasi-steady states are large in comparison with the observing 
time. Here we present a new set of stationary wind models and obtain their associated X-ray luminosity, surface brightness and spectroscopical 
temperature. Certainly, the presented X-ray models do not include short-scale variations associated with interacting binaries. These are 
associated with the point-like X-ray sources, whose luminosity must be removed from the diffuse X-ray component (Moffat et al. 2002). 

Consider the Chevalier & Clegg (1985, CC85 henceforth) fast superwind model, which assumes that the flow is steady, spherically- 
symmetric, adiabatic and ideal. Each star cluster can be defined by a set II* = {R S c, q e , {<7m, Vooa}} containing 3 parameters: the cluster 
radius, R sc ; the energy deposition rate per unit volume, q c = E/V sc ; and the mass deposition rate per unit volume, g m = M /V sc - Alterna- 
tively, the adiabatic terminal speed Kx>a = (2q e /q ln ) 1 ^ 2 may be used. For the interior of the cluster, r < R sc (region A), the equations of 
conservation of mass and momentum are 

(1) 

and 

du dP 

pu— = - q m u, (2) 

dr dr 

respectively. Above, only models with finite central densities (p c > 0) are considered, i.e. it is a requirement that u c = u(0) = km s -1 . 
For region B, the corresponding equations are 

= M. (3) 

4nur 2 

and 

du dP 

pu— = - — . (4) 

dr dr 

For both regions, the conservation of energy can be stated as a steady state Bernoulli-like equation 

e=\u 2 + ^-- = \vl A , (5) 
2 7 — 1 p 2 

where e is the total energy per unit mass. Since the gas is adiabatic and perfect, 7 = 5/3. Let's consider now a polytropic perfect gas with 
equation of state 

P = Kp^ = nkT, (6) 

where 77 is the polytropic index and K > is a proportionality constant. This equation replaces l(5} in the polytropic case, where the entropy 
variations are parametrized through a constant 77. 

The combination of I0 and (|6) yields the explicit form of the energy conservation per unit mass for region B 

1 2 p 

£bp = — u + (77 + 1) — = constant. (7) 

2 p 

As expected, the polytropic and the CC85 solutions have the same algebraic structure on r ) _R SC . They are isomorphic through the 
transformations 7 — > (77 + l)/»7 and ebp = (1/2)V^ P (V00A — * Voov). The constant Voop is arbitrary and is related to the temperature 
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Figure 2. Panel (a): comparison of the adiabatic isomorphism velocity profile for region A, equation {§}, with that of the pseudo-adiabatic solution, equation 
{16}. In general, the solutions are alike for the same set of parameters. Here, r) = 3/2 and Voop = VqoA- Panels (b), (c) and (d): normalised velocity, density 
and temperature profiles for pseudo-adiabatic polytropic solutions with r\ = 3/2 (solid lines), r\ = 1 (dashed lines) and rj = 2 (dotted lines), respectively. 
The temperature profile is by far the most sensitive to changes of r\. For superwinds with terminal speeds of the order of Vg 1, T sc a > 10 7 K. 



at r = R sc , or, equivalently, to the speed at which the superwind leaves the star cluster surface. However, this arbitrariness vanishes when 
additional constraints or boundary conditions related to the parameterized effects are imposed on the model, for example, VooP < VooA when 
just dissipative processes are involved. When more specific information about the nature of the processes is available, explicit expressions 
for the constraints can be given as threshold lines. For region A, the isomorphism is not exactly the case; nevertheless, as discussed below, 
equation ((7} remains a good approximation there. 

A continuous gas acceleration requires of the following condition at r = R BC 



(V + 1) Ps< 



-,1/2 



< Vex 



The isomorphic solutions are self-similar if written in terms of the dimensionless variables R = r/R BC and U = u 2 /V^ P 
R = D 1 U 1/2 [l + (6r ] + 5)U]-£$ , r<R sc , 
R = D 2 U' 1/i (l-U)-^ , r^7? sc . 
Accordingly, the boundary condition l[8j becomes 

U sc = 



1 



277 + 1 

From equations $9$ and < | 1 Ob and condition dl It . one can find that for an accelerating solution 



£>i = (2*7 + 1) 1/2 
and 



(2r?- 



-1/4 



1 + 



6*7 + 5 
277 + 1 



<>il-r. 



77/2 



27,- 



(8) 

(9) 
(10) 

(11) 
(12) 
(13) 



The actual polytropic solution for the innermost region (r < R sc ) can be obtained from one of the pfaffian forms related to the 
conservation laws. The simplest pfaffian form is 
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(Ap 2+ v - r 2 ) dp + 4pr dr = 0, (14) 

with A = 9(r)+l)K/ (qmV)- The associated ordinary differential equation is inexact, but it has a unique integrating factor: If = p'. After 
multiplying by it, the next solution can be found with ease 

r 2 +^ P ^+C P ^ = 0, (15) 
3r] + 2 

where C is an arbitrary constant. Such a solution has several branches according to the values that C, K and 77 take. By using equation 10 and 
the transformation epp = f V^p as closure relations for the the true polytropic solution given by equation J 1 5b . a simple parameterization 
that we call the pseudo-adiabatic solution can be obtained 

# = Dip(y 1/2 (l + 3(277 + l)[fr 2/3 (r<R EC ). (16) 

In what follows, we shall consider the polytropic solutions l !16t and J 1 0| > satisfying the boundary condition i ll It . Thus, 
Dip = 2 4/3 (2t7 + 1) 1/2 (17) 
and the values of the hydrodynamical variables at r = _R SC are 

p sc = (2r ? + l) 1/2 ^i, (19) 



3V« 



and 



where pi is the mean mass per particle for a completely ionized interstellar plasma and k is the Boltzmann constant. For solar abundance 
Pi = 14/23pH- The central values are 



1 Pi 



T * = WTT)T V ^> (21) 

Pc = Dip— — , (22) 

OV00P 

and 

Pc = D 1P ( ^SlV^. (23) 



k 2»7 + l. 

Notice that the models are valid just for 77 > 0. For simplicity and in order to cover the polytropic states of astrophysical interest, we 
shall further consider that 77 > 3/10 in region A. The isomorphism approaches the pseudo-adiabatic solution for large 77 (the isothermal case). 
For intermediate values of 77 they are similar [see panel (a) of Fig. [2] since they predict the same values for the hydrodynamical variables 
at R sc and the same central temperature and velocity. Their densities and pressures just differ by factors of D± / Dip at the centre and their 
ratios reduce towards i? sc . 

Since the effective polytropic index depends on the balance of heating and cooling among other factors, our models can parameterise 
these effects when complemented with acceptable boundary conditions or constraints, as we discuss next. In principle, zones with distinct 
entropy variations and properties can be modeled using a piecewise polytropic index, or, in the limit, a slow varying one. For the case of a 
continuous accelerating solution, equation implies that any change in the available energy per unit mass is concomitant with a change of 77 
in the same direction. Nevertheless, notice that velocity discontinuities produced by shock fronts are not prohibitive. For a smooth accelerating 
solution with parameters II, , models with 77 < 3/2 have higher velocity profiles and lower density profiles than the corresponding adiabatic 
solution (77 = 3/2, VooA = Voop)- On the other hand, their temperature profiles have a bimodal behaviour: in region A they are higher, 
whereas in region B they are lower and very sensitive to the adopted value of 77. Models with 77 > 3/2 have the opposite behavioufl (see Fig. 
|2j. In the external zone, all the density profiles fall off asymptotically as M / A-kVoopt 2 . In connection with the radiative case, Silich, Tenorio- 
Tagle & Rodriguez-Gonzalez (2004) have shown that whilst the velocity and density profiles are barely affected by the small variations of 
the effective terminal speed promoted by radiative cooling, the temperature profile can be changed drastically at distances larger than a few 
times the star cluster radius. Given that for dissipative solutions V00P < VooA and that the relative behaviour of the profiles has to scale 
accordingly, it follows that models with 77 > 3/2 are appropriate for region A, whereas those with 77 < 3/2 may be used for region B. In 
Fig. [3] we reproduce the numeric results of Silich et al. (2004) for radiative winds using the aforesaid method. At least 3 different ?7's were 
required for each case presented. The profiles were constructed using the X-ray luminosity (Section l2.lt to estimate the energy radiated by 



1 Notice however that a very large rj produces very slow winds which are susceptible to instabilities and might have different driving mechanisms. 
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Figure 3. Silich et al. (2004, fig. 4) radiative temperature profiles for clusters with II, = {R sc , P c = 10, E 38 = {500, 2000, 3000}, V 8 = 1}. They were 
here constructed using piecewise continuous polytropic laws. For the cluster cores (region A), 77 ^ 3/2 for the three mechanical luminosities. The respective 
free winds have different thermodynamical configurations: Bo represents a quasi-adiabatic wind with E-jg = 500, 771 £S 1.45 (dashed-line) and r?2 ~ 1.15 
(dash-dotted line); likewise, B\ is a radiative wind with Ezs = 2000, 771 ~ 1.3 and 772 ~ 0.42; and B2 is a strongly radiative wind with Ess = 3000, 
771 ?3 1.19 and 772 ~ 0.22. In the later case, the velocity-radius relationship adopts 77 = 3/2 for T > 10 4 K (dotted line). The assumed metallicity is Z = 
Zq . The squares show the external temperature profile expected from the adiabatic solution, which better applies for clusters with smaller masses (mechanical 
luminosities). 



the cluster core, since Ax ~ A tota i for temperatures in the order of ~ 10 7 K and Z ~Zq. The emissions of the external sub-zones were 
weighed (see Appendix [At using the total cooling function tabulated by Plewa (1995). From this, the respective changes in K»p were found, 
and with the aid of 10, also the corresponding values of r\. As a conclusion, it can be pointed out that just a change of r\ (adequate for each 
zone) is sufficient to reproduce and explain the radiative profiles. For region A, we found that 77 is always close to 3/2. However, a closed 
analytical estimation to find explicitly the adequate values of 77 for the free wind is still lacking. In principle, just a knowledge of the star 
cluster parameters and the overall cooling function is sufficient to obtain it. A general discussion of the radiative case will be presented on a 
forthcoming paper. 



2.1 X-ray luminosity 

The cumulative X-ray luminosity of a spherically symmetric gas distribution is given by the integral 

L x (r)=4n [ n(r') 2 A x [T(r'),Z]r' 2 dr', (24) 



where r is the radial coordinate from the object centre, r' the corresponding dummy integration variable, n is the particle number density, 
and Ax is the X-ray emissivity function, which in general depends on the plasma temperature and metallicity. Here, the later is assumed 
homogeneous; nevertheless, tools to deal with metallicity gradients are given. 

The actual dependence of the superwind X-ray emission on the star cluster parameters is determined by the functional relation of Ax 
with the hydrodynamical variables. This functional dependence synthesizes the physics behind the emission mechanisms. In Appendix [AT] 
we give an initial approximation that corresponds to a constant emissivity, whereas in Appendix lA2l we give transformation laws that are used 
to show how to overcome such a limitation. Those laws can also be used to incorporate the effect of metallicity gradients or non-collisional 
ionization equilibrium. Below, we concentrate on region A because, as is demonstrated in Appendix I All region B contributes less than 20 
per cent to the superwind total X-ray luminosity. 

In Appendix IA2I it is demonstrated that for a realistic emissivity given as piecewise continuous power laws, Ax = A a T a , the 
cumulative X-ray luminosity is given by 

L XA ( V ,{a},U) = X A U s/2 Y / ^Av,a,U)A a T c °', (25) 

a 

where 

Aa = — -577^- p 1 ( 26 ) 
350 

F A (r,,a,U)= I i v+ ) -£/)*, (27) 

My, a) = 70(a + 2877 + 14), (28) 



On the X-ray emission from massive star clusters and their evolving superbubbles II 1 




Figure 4. Panel (a): region A cumulative X-ray luminosity normalised by a factor of X\ /Df p , which contains the dependence on the cluster parameters. The 
values of the polytropic index are r\ = 3/2 (solid line), r\ = 1 (dashed line) and r\ = 2 (dotted line). Panel (b): ratio of the X-ray luminosity of metals (solid 
lines, region A) to that of hydrogen for different core temperatures (Tkev). Here r\ = 3/2. The dashed line is used as a reference. From bottom to top of the 
panel, the metallicities are 0.1,0.2,0.5,1.0,2.0 and 5.0 Zq. The axes marks are logarithmically spaced for a better visualization. 



and 



<j> 2 (r), a) = [17a 2 + a(501 + 952??) + 9528^(1 + n) + 2382]. 



(29) 



In these expressions a G [—4, 4] and U £ [0,1/(2^ + 1)]. 

To account for a realistic emission, we assume collisional ionization equilibrium (CIE). This assumption is valid for the central region, 
where the ionization and recombination time-scales used to assess the collisional equilibrium (Mewe 1984, 1997) are smaller than the 
radiative and adiabatic expansion cooling time-scales (see Ji, Wang & Kwan 2006). Outside of the cluster, CIE does not hold for r > 2-3i? sc 
(T < 3 x 10 6 K); however, the predicted X -ray emission of such region is negligible. Hence, we use the Strickland & Stevens (2000, SS00 
hereafter) X-ray emissivity tables, which are based on a recent version of the Raymond & Smith (1977) code. These tables give the X-ray 
emissivity as a function of temperature for two energy bands appropriate for Chandra and XMM-Newton observations: a soft band that goes 
from 0.3-2.0 keV and a hard band from 2.0-8.0 keV. They account for continuum emission by incorporating the free-free (dominant for T > 
3 keV), the free-bound and the two photon (significant for T ^ 2 keV in the energy range 1-10 keV) processes. They also account for the 
line emission, which dominates the spectrum at temperatures around and below 1 keV. Additionally, the contributions from hydrogen and 
metals were separated to allow the inclusion of arbitrary metallicities. Therefore, the total X-ray emissivity can be written as 



A X (T, Z) = Axh(T) + ZAxm(T). 



(30) 



We fitted power laws and rational functions to the SS00 tables using a trust region method (see Table lD2l and Table lD3l . As an outcome, 
it was found that for fast winds ( Vooi> > 1000 km s _1 and T BC > 9 x 10 6 K) the hydrogen component in both bands well resembles a 
bremsstrahlung law, whereas the emissivity of metals has an almost constant value at high temperatures (2.6-8.6 keV) and a sharp power 
law behaviour at the other end (0.8-2.6 keV). Using these results (marked with ** in Table lD2l and evaluating formulas i25\ and l !27t at 

U S c = 1/(277 + 1) we g et that 



ixA.total = XaU sc 



3/2 



^(r?,-,t/ S c)A H T c 1/2 + Z 



FAv, -5, E4c)AmiT c - 5/2 + T A (j], 0, U S c)Am 2 T° 



(31) 



Bearing in mind equations J2 1 b and 

LxA.total = ( 5 - 12 x lO^ergs- 1 ) 
or, in terms of Tk c v, 

LxA.totai = (9.31 x lO^ergs" 1 ) 



£■38 



and taking 77 = 3/2 we have that 



1 ^sc,pc ' g 



-flsc,pc-' kcV 



1.00 + 1.32Z 



2.05 

T 3 

-'keV 



1.47 



1/2 



(32) 



(33) 



where the first and second terms between square brackets give the relative contribution from hydrogen and metals, respectively. Since a term 
oc T 1 / 2 was factorized, the bracketed quantities show the variations of the SS00 emissivity function with respect to bremsstrahlung. The 
X-ray luminosity of warmer (T sc < 9 x 10 6 K) and generally slower winds (Vs < 1) can be obtained in a similar way using equations 
l|25]l-l[29]l together with Table |P2l 
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In Paper I, the total X-ray luminosity of regions A and B was estimated by approximating the temperature and density by the constant 
values that they take at R sc according to the CC85 adiabatic model, i.e. by setting p = p sc and T = T sc with rj = 3/2. Consequently, it was 
found that LxAB.totai = 47rapAx(T sc , Z)p 2 c R^ c /3pl. Explicitly, this is 

Lx.ab = (1.27 x 10 34 erg s" 1 ) A ^( T ^ Z J E ™ = ( L08 x if/Vg s^ 1 ) ^,^3(^)^8 (34) 

In the above formula, the effects of the actual thermodynamical structure of zones A and B and the realistic X-ray emissivity (SSOO) were 
incorporated through the fiducial coefficient a p ~ 2, which smoothed the differences between the analytical and the numerical calculations 
that were carried out for winds with Vs ~ 1. Because such calculations also showed that the contribution of region B was always less than 
one-quarter of that of region A, the former was also incorporated into d34t through a p . These results agree fairly well with equations ( IA6b 
and JA12t which just assume a constant X-ray emissivity. In contrast, equations {32} and J33 b explicitly show their physical dependence 
on the thermodynamical structure and the SSOO X-ray emissivity function through the numerical coefficients obtained, the separation of 
the contributions of hydrogen and metals, and their functional dependence on the terminal speed and temperature, respectively. Another 
advantage of these formulae over that of Paper I is that they can be used for very fast winds with V00P ~ 2500-^1500 km s - 1 . 

Additionally, in the current formalism, equations J25t and d 1 6b allow to calculate the cumulative X-ray luminosity profile of region A for 
different polytropic indexes. If it were needed, the luminosity of region B could be calculated using our transformation laws. The cumulative 
X-ray luminosity profile of the cluster core for different values of 77 is presented in panel (a) of Fig. [4] In the same figure, panel (b) shows 
that for Ikcv ^ 3 keV the contribution from hydrogen is comparable (within a factor of ~ 2) to that from the metals when Z ~ Z©, whereas 
for lower metallicities hydrogen can dominate in the same temperature range. On the other hand, metals dominate the cluster core emission 
at lower temperatures when Z ~ 0.5-1 Zq or higher. In general, hydrogen always dominates for Z < 0.2 Z© and metals for Z > 2 Zq. 
We remark that metallicity gradients or a different X-ray emissivity can be handled using the aforementioned transformation laws and the 
adequate polytropic index. For instance, when the CIE assumption does not hold one can use non-equilibrium emissivity functions like the 
one presented by Gnat & Sternberg (2007), whom for solar metallicity found that A inoncq oc 5.6 x io- 20 T-°' 46 erg s 1 cm . When CIE 
does not hold inside of the cluster, formula J25b can handle this new emissivity function with ease. By separating components and because of 
their ability to handle arbitrary emissivities, our formulas round out the X-ray luminosity model for the gas contained within the star cluster. 



2.2 Surface brightness 

To calculate the surface brightness and other projected quantities we use the Abel transform (see, e.g. Yoshikawa & Suto 1999) 



Q(s) = 2 q(r)-= =^dr, (35) 

is V ^ ^ 

where r has the same meaning as before, s = r^B is the projected radius, r<j is the angular diameter distance to the object and 6 the projected 
subtended angle. The factor of 2 takes advantage of the spherical symmetry (see Fig. [5}. Since we are dealing with finitely localized gas 
distributions, the traditional upper limit (00) is replaced by the size of the region of interest. This size is determined by the gas temperature. 
We consider an X-ray cut-off temperature of T cut = 5 x 10 5 K for all regions. The corresponding radius is R cu t- 

The surface brightness can be obtained by setting q = n 2 Ax in d35t . For the superwind, it is given in terms of elliptical integrals in 
non-canonical form, the reduction of which requires the awkward calculation of the roots of several quartics for each s. None the less, the 
whole superwind (region A+B) surface brightness can be fitted by a function of the form 



ctab(s) = Aae 



1 



(36) 



For r\ = 3/2 we have that «i w 3 and Q2 ~ —5/2. The associated central brightness can be found in terms of hypergeometric functions. 
For fast winds, Vs > 1, it is given by 



A x (T ac , Z)E 2 34 , 2 Ax,-23(T BC ,Z)E; 



^ab - 1.28 - — — = (1.9 x 10- ergs- (37) 

/- i n- f * J sc'oo - f1 'sc,pc- i kcV 

From this formula and Fig. [6] it becomes clear that even in the adiabatic case (rj — 3/2) the superwind X-ray surface brightness drops 
extremely fast (oc s _15//2 ) for s 3> R ac - Thus, although in principle clusters in the quasi-adiabatic regime per se can generate very extended 
X-ray emitting haloes, see equation l lA7t . their surface brightnesses far away from the cores could be too weak to be detectable, specially 
when the integration times are short (see Lang et al. 2005). 



2.3 Weighted temperatures 

Several definitions of the projected X-ray temperature exist, that is, of the temperature Tx by which the spectrum of an often multi- 
temperature plasma is fitted as if it would have just one component. The most frequently used ones are of the form (see e.g. Navarro, 
Frenk & White 1995) 
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Figure 5. Abel transform of a spherically symmetric function with compact support on a ball of radius R and projected angular diameter 6. For r d ^> R the 
Radon Transform lines (dashed), approach the Abel Transform lines (dotted) since Sfi = r d sin(#) r d and s_4 = r d tan(S) rs r d 0. The right angles 
that the projected radii and the respective projection lines subtend are indicated with black and white squares. The spherical symmetry of the function warrants 
that the integral yields the same value along both lines of sight. The difference is that in the astrophysical case the distribution is projected along the 1Z lines, 
i.e. the values at points p 1 are mapped to points p. Compare the difference in the auxiliary projection plane. 
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Figure 6. Normalised X-ray surface brightness of a superwind (r/ = 3/2). Soft line: actual surface brightness, dashed line: fitted profile. Note that although 
the outer layers of the core contribute relatively more to the central X-ray luminosity (Fig. [5), no limb-brightening is produced there. The overall superwind 
surface brightness sharply drops in the external zone. 
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f n 2 ATdV 

Tx = W' (38) 

where A is a weighting function that may be the X-ray emissivity function (see Mazzotta et al. 2004 and references therein). The emission- 
measure weighting scheme uses A = 1, whereas the so called (bolometric) emission weighting scheme assumes that A = A(T) oc T 1 / 2 , i.e. 
it assumes just free-free emission, notwithstanding the fact that this later mechanism dominates just for T > 3 keV. Mazzotta et al. (2004) 
have shown that theoretically, none of these definitions is an accurate approximation of the spectroscopic temperature of multi-temperature 
plasmas. However, they have argued that since noise, instruments response, and instrumental and cosmic backgrounds among other factors 
affect and distort the observed spectrum, sometimes it is possible to fit a single-temperature plasma model to the data. In their work, they 
simulated the expected spectra of galaxy clusters (Z =Z©) taking into account the sensitivity of the Chandra and XMM-Newton instruments 
in the soft band. They found that for T ^ 3 keV the synthetic spectra (with temperature T spcc ) could be reproduced -at a level better than 5 
per cent- by weighting with a function of the form A = A(T) oc y -3 / 4 . Such a function yields a spectroscopic-like temperature Tsl that 
can be used as an estimator of the actual temperatures of the isothermal models that might be fitted to real spectral data. 

None of the above schemes take into consideration arbitrary metal abundances. Vikhlinin (2006) sought to extend the previous result for 
T > 0.5 keV and arbitrary metallicities. Although his algorithm yields a significant accuracy gain when lower temperatures are involved, it 
does not provide analytical estimations because his weighting function has to be constructed ad hoc (it depends on the measurement device 
characteristics and individual observational conditions). We present here a theoretical weighting scheme based on the SS00 tables. Since 
we do not consider observational artefacts, our spectroscopical temperature Tx corresponds to the ideal projected X-ray temperature of the 
models. Nevertheless, our scheme can incorporate artefacts easily, provided that the overall weighting function could be parameterized in the 
same way as the one found by Mazzotta el al. (2004), i.e. as a power law (or series) in T. 

For bubbles and superbubbles we give preference to the projected temperature map as a diagnostic tool 

_ Jn 2 AT(r 2 ^s 2 r^ 2 rdr 

Tx(s) ~~ fn 2 A(r 2 - S 2)-i/2 rdr • (39) 

For superwinds, the emission-measure and the Chandra and XMM-Newton spectroscopic-like temperatures (whether the later applies or 
not depends on r\; for r\ « |, it is a valid measure when Vs > 1.65) can be obtained from equations d25b and d27b 

T X i(r/,a, Use) = — F -, TT z = — ■=—, TT > T c , (40) 

Lxa(?7,«, Use) •>A\J)t U 3 c) 

where a — for the first estimator and a = —3/4 for the second one. Finally, our X-ray emission weighted temperature is 

lx2(r],a, U S e, Z) — — — — , (41) 

LxA{ri,{a\,U S e) 

which depends implicitly on metallicity. For superwinds with Vs > 1, rj > 3/2 and Z =Z©, all the above estimators yield that Tx ~ 0.95 T c 
and that Tx gets closer to T c as r\ increases. There is just a slight difference at the other end: for r\ w 0.3 we obtained that Tx ~ 0.85 T c . 
This implies that the cluster central temperature is adequate for representing the spectroscopical one. For superwinds, T sc /T c G [3/8, 1), so 
in many cases, the temperature difference between the centre and the edge might not be significative. However, for superbubbles -as we shall 
see- such a difference does matter. 



3 BUBBLES WITH SHELL EVAPORATION 

On the theoretical side, X-ray models for superbubbles with evaporation from the cold, dense and thin outer shell (W77, region C in Fig. 1) 
have an interesting analytical behaviour. During the sweep-up phase, thermal conduction fluxes across the contact discontinuity promote the 
injection of radiatively cooled material into the cavity that contains the hotter (shock-thermalised) but much more tenuous superwind. As a 
result, the gas distribution has very steep density and temperature gradients and to account for the actual shape of the emissivity function 
becomes crucial. On the other hand, in some cases the observations pose another challenge, for after comparing with the C95 standard X-ray 
luminosity model, discrepancies of up to two orders of magnitude have been found (Chu, Gruendl & Guerrero 2003; Dunne et al. 2003). 

Most analytical models for spherical bubbles on the sweep-up phase (Castor et al. 1975; W77; Hanami & Sakashita 1987; Koo & McKee 
1992a,b; Garcfa-Segura & Mac Low 1995) share the following similarities: 

(i) In all models, the compact outer shell is driven by the thermal pressure of region C. There, the shocked wind is isobaric (except in a 
small zone near the secondary shock) and its pressure equals that of the shocked ISM. Therefore, the pressure of the hot bubble is only a 
function of time. 

(ii) Similarly, it is demonstrated that the bubbles have a pseudo-adiabatic interior because the characteristic cooling time of the shocked 
superwind is larger than the bubble dynamical time scale (see, e.g., Mac Low & McCray 1988 and Koo & McKee 1992a,b). The prospective 
evaporation of cold gas from the external shell does not modify this behaviour significantly, even after the evaporated mass dominates the 
region. 
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(iii) Under this common premises, but considering different additional physics, the authors determined the bubble geometrical evolution 
and its energy budget, finding explicit expressions for the principal and secondary shock positions, the pressure, and the total energy: R p (t), 
R s (t), Pc(t) and E c (t), respectively. 

(iv) With exception of the Koo & McKee model, which just considers the intrinsic self-similar evolution of region C, the rest of the 
models assume that the thermodynamical structure is determined by the thermal evaporation of cold material from the external shell. The 
later process is governed by the classical thermal conduction theory (Cowie & McKee, 1977). The predicted temperature profile is of the 
form 

/ s 2/5 

T c (r) = T Cc y--^-) ■ (42) 

(v) The density profile is derived from the isobaric condition stated in (i) 

-2/5 

Rr, 



nc(r)=ncc(l--^-l , (43) 



where Tec and ncc are the extrapolated central temperature and density. Note that the above profiles are spatially self-similar. 

The points above imply that the dynamical evolution (linked to the motion energetics) and the thermodynamical structure are totally 
independent, regardless of the choice of the bubble model. Each model assumes different physical conditions that lead to different evolutions: 
W77 considered constant energy injection rates and ISM densities. They showed that if radiative cooling (when the bubble interior departs 
from a quasi-adiabatic state) and distinct ambient pressures are incorporated into their model, the bubble size can change but the self- 
similarity of the thermodynamical profiles is not affected. Hanami & Sakashita (1987) and Garcia Segura & Mac Low (1995) assumed 
different ambient stratifications: they considered ISM densities that fall off as po oc and p oc r~ 2 , respectively. The later authors also 
incorporated variable mechanical luminosities. Again, the geometrical expansion is different, but the shapes of the density and temperature 
profiles are not disturbed. 

The coupling of the physics resulting from (iii) and (iv) is done through Tc c and nc c which inherit the time dependence of the geo- 
metrical model, i.e. they depend on the derived R p (t), R s (t), Pc(t) and Ec(i). They also enclose the efficiency of the evaporative process 
given by the numerical value of the thermal conduction coefficient (Spitzer 1956). It can be shown that as long as conditions (iv) and (v) 
are satisfied, the X-ray emission and other integral properties like the evaporated mass (which dominates in this zone) are independent from 
the geometrical (dynamical) model. These considerations are taken into account in the X-ray models here presented, where they are used 
to derive properties that are valid for all evaporation dominated models. Such properties can be helpful to constraint the related physics and 
serve as auxiliary tools at the time of analysing the observations. 

For simplicity and in order to explore the parameter space of a specific model, we shall consider that of W77, unless otherwise stated. It 
that model, the leading and reverse shock positions and the central temperature and density are given by 

R p = (67pc)£& l V /B *r. (44) 

/ F 3 \ 1/10 

* = (26 P c)Uf!) #*> ^ 
T Cc = (5.25 x 10 6 K)£f 8 /35 n 2 / 3B t 6 ~ 6/3B , (46) 
and 

n Cc = (1.63 x 1CT 2 cm- 3 )4 6 8 /35 nJ 9/3 V 2/35 , (47) 

where no is the particle number density of the ISM. For normal OB associations, McCray & Kafatos (1987) have pointed out that initially 
(at ages up to 5 x 10 6 yr), the injection of energy (a few times ~ 6 x 10 erg s _1 ) is dominated by the most massive but much less numerous 
O type stars which blow bubbles with sizes ~ 100 pc. At later times and up to 5 x 10 7 yr, the injection rates are sustained by the supernovae 
explosions of B type stars; specially, after the most massive star have died. The supernovae inject energy at a steady rate of N*a x 10 35 erg 
s _1 (where a is a constant of order unity and iV* is the total number of stars with masses above 7 Mq) and drive a new outflow over the 
previously injected gas. For rich OB associations, the wind dominated stage can be as important as the supernovae dominated one, i.e. for 
associations or clusters with E-^s > 1 (Abbott, Bieging & Curchwell 1981). When these effect are enhanced, models like those of Hanami & 
Sakashita (1987) and Garcfa-Segura & Mac Low (1995) might be more adequate to describe the evolution of the corresponding superbubbles. 
Here however, we follow McCray & Kafatos (1987) and assume that the energy input rate remains constant through time in order to analyse 
later the case of M 17. 

3.1 X-ray luminosity 

In its simplest form, the total X-ray luminosity of a well-developed bubble with shell evaporation is 

ixctotai = 47rn 2 3c i? 3 A x (^cut, Z c ;T Cc ), (48) 
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Figure 7. Solid lines: evaporatively averaged soft band emissivities resulting from the C95 constant value and the SSOO tables (hydrogen: Axh, metals: Axm)- 
Here, T7 is Tq c in units of 10 7 K and Zq =Zq. The lines for the C95 model and Axm scale linearly with metallicity. The former was extrapolated to the 
temperature range 10 6 -2.5 X 10* K. The best-fitting approximations to Axh an d Axm are shown as empty squares. Significant deviations from the C95 
predicted values are observed as Tq c 10 6 K. The deviations induced by the reverse shock when R s /R p = 0.6, 0.7 and 0.8 are shown as dashed-dotted 
lines cascading down from the respective averaged emissivities. The axes marks are logarithmically spaced. 

where 

A x = | f tf 1/2 (l - tf 5/2 ) 2 A x (tf, Z C ; T Cc )M, (49) 

■3 = Tc/Tcc is the normalised temperature and i? cut = Tcut/Tbc is the normalised cut-off temperature. Generally, Ax (on the inte- 
grand) cannot be put in dimensionless form with respect to ■& and, as a consequence, it will depend on Tec as a parameter. The variable 
Ax(#cut, Zc\Tcc) can be thought as an "evaporatively averaged" X-ray emissivity. It is practically independent of the geometrical evolution 
when R s /Rp < 0.4 (see Fig.|7J. This is what we call a well-developed bubble and thus the upper limit of the integral in J49l > is justified. For 
0.4 < R s /R p 0.6 there is just a slight geometrical dependence for bubbles with Tq c < 3 x 10 7 K. If this ratio increases so does the later 
effect, specially in the ~ 10 6 K temperature range. Because in all current models bubbles get well-developed at early times, we shall just 
discuss that case. 

In Paper I, the bubble X-ray luminosity was estimated using the C95 model which assumes that for the Einstein satellite Ax(T, Z) — 
ZqAo = 3Z x 10 -23 erg cm 3 s _1 on the temperature range 10 6 -10 7 K. Therefore, one obtains that Ax,C95 = ZcAo/($ C ut) with 
J(tf) = (125/33) - 5tf 1/2 + (5/3)i? 3 - (5/ll)i9 11/2 . However, in some cases, the differences that arise when one considers the more 
realistic SSOO emissivity function can be substantial, as shown in Fig. [7] There, the average of the SSOO emissivity given by equation ( 1301 ) 
displays separated the expected soft-band X-ray contributions from hydrogen and metals. The corresponding X-ray luminosity is 

ixctotai = 4-Kn 2 Cc R 3 p [Axh(Tcc) + ZcAxm(Tcc)] . (50) 

The above formula is just corrected for the effect of the shape of Ax; the variables that depend on the geometrical evolution (R p , Tec 
and ticc) are still undetermined as functions of the problem parameters (no, E, Voop) and time. Such functional relation is a consequence 
of the assumed bubble model. However, a great advantage of ( 150b is that it just depends on variables that are either directly observable or 
spectroscopically derivable. It is clear that the key issue is to find relationships between the central values of temperature and density and 
their spectroscopic estimators (Section[33J. 

On the purely theoretical side, if one also assumes the geometrical growth given by W77, an analogue to the C95 X-ray formula follows 

Lxc,W77 = (3.25 x 10 35 erg s- 1 )AxEH /35 n 1 7/35 t 1 6 9/35 , (51) 

where Ax is the average of the SSOO emissivity. Since it is smaller (for all T) than the constant value assumed by C95, the bubbles evolve 
even more adiabatically and the use of the W77 model is a safe assumption. This ensures complete consistency with the emissivity function 
that we are using and our aim of incorporating additional physics. 

3.2 Surface brightness 

For a single power-law X-ray emissivity, Ax = A a (Zc)T a , the surface brightness was decomposed (Appendix|Bj into three components: 
a temporal component related to an intensity amplitude, Ac(t); a weighting emissivity that implicitly varies with time because of its depen- 
dence on Tec, A Q = A aj _23Tc c (though \ a = Ao = 3Zc is a constant for the C95 model); and a spatial profile that slowly varies with time, 
Sc(s) [or as a function of the normalised temperature]. In terms of these components ac(s,t) — Ac(t)X a Sc(s). Their respective 
model-independent expressions are given by equations dB 1 t - dB3t . 
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For the W77 evolutionary model the amplitude factor is given by 

A { i nc . in3! -1 -2 \ £il9/35 31/35,-23/35 , cos 

■Ac(t) = (1.05 X 10 ergs pc )E 3S ' n ' t 6 ' . (52) 

Because of the thermodynamical structure self-similarity, the shape of the spatial profile (Sc) in terms of the bubble normalised projected 
radius (s/R p ) depends solely on the the instantaneous value of the central temperature, not on the dynamical evolution. This implies that for 
a fixed Tc c such a shape is the same for all bubble models accounting for shell evaporation. Note that A Q is also fixed for a given Tc c and 
that the prescription of a particular bubble model merely implies a scaling of the intensity amplitude. Thus, the results here presented are 
very general and can be applied to the overall time evolution of a single object as well as to different objects with distinct parameters at a 
fixed age. As we shall see, the form of the assumed X-ray emissivity will play an important role in the subsequent discussion. 

Approximations for 5c as a function of temperature were implemented in equations dC9t and JC lit for a single power-law emissivity 
and the realistic SS00 emissivity, respectively. The first one allows a straightforward evaluation, while the last one, more elaborated, is 
necessary for a more robust discussion. In terms of s £ [0, R C ut] the former becomes 

I SlT l + 2a -i ff=max(-p— 

ScA{s,a;^,iO=[l-—j [(l -£)-!- ^(l-M-fca)]^^ P P , (53) 

where the square-bracketed quantity written in terms of the dummy variable £ is meant to be evaluated at the indicated limits (as the result 
of a definite integral would), A indicates that we used a pseudo-triangle approximation (Appendix let. max(s / R p , R B /R P ) incorporates the 
reverse shock position and J~c,c (t, Tf ; a) given by equation JC8t is an algebraic expression that involves the arcsin function. The central value 
of the spatial profile is obtained by evaluating the above equation at s = 0. The upper limit for £ then becomes max(0, Rs/R p ) = Rs/R p 
and we have that 



iSca(0, a; 



l + 2a 



Tcc) \Tcc 



(54) 



where Tr s = Tc c (l — Rs/R p ) 2 ^ 5 is the temperature at the reverse shock position according to some specific bubble model. When _R a is 
ignored (which is equivalent to putting i? s = 0) as in the C95 X-ray luminosity model, Tr. s /Tc c in the previous formula becomes unity. 

As a first step, we obtain the surface brightness corresponding to the C95 X-ray luminosity model. Theoretically, this model is equivalent 
to considering that the X-ray emission is proportional to the emission measure, i.e that Lxc = ZcAqEMc (a = 0, Ax = ZcAo). 
According to this, one has that while interstellar bubbles present a practically flat profile, superbubbles driven by the most massive SSCs 
exhibit a noticeable limb brightening near the contact discontinuity due to their more powerful mass and energy injection rates, which 
translate into higher central temperatures. This is shown in panel (a) of Fig. [8] The same panel also shows the direction of the spatial profile 
expected time evolution for a fixed set of model parameters. As time progresses, the profile becomes flatter and weaker due to a continuous 
and slovjf] decrease of the central temperature. However, the energy injection rate is not fixed as it increases by a factor of 3 after 3-Myr 
(Leitherer & Heckman 1995). This implies that a slight limb-brightening increment is expected for both bubbles and star clusters at this 
age, when the first supernovae start to explode. These results derived from our analytics and the ones presented below (Fig. [5} can also be 
used to reconstruct the simulated X-ray brightness profiles of Wrigge et al. (2005) for the Wolf-Rayet bubble NGC 6888. In their work, the 
Garcfa-Segura & Mac Low (1995) bubble model was used and the required X-ray properties were obtained by ad hoc adjusting the thermal 
conductivity coefficient in order to vary Tc c - We remark that the different stratification of the ISM that they considered does not produce a 
divergence from the analytical shape of our spatial profiles, since as it was previously stated, the thermodynamical configuration of region 
C is essentially the same. There are however two conspicuous differences with respect to the W77 model: in their case a) the spatial profile 
evolves faster (Tcc oc t~ 2//7 and nc c oc t~ 12 ^ 7 ) and b) the predicted size of the bubble is not the same. 

In panel (b) of Fig. [8] the C95 model is compared with one that assumes an emissivity of the form Ax = ZqAcT 1 ^ 2 , i.e. a 
bremsstrahlung dominated emission, which is a commonly used approximation. The main difference between them is that the second one 
practically does not show any increment in surface brightness near the contact discontinuity. This is in better agreement with observations of 
the superbubbles associated to the Rosette and Omega Nebulae (Chu, Gruendl & Guerrero 2003; Townsley et al. 2001, 2003) for which no 
limb-brightening was detected. One then might suspect that the sharp peak predicted from a constant emissivity could be an artefact. If an 
emissivity that decreases towards the cut-off temperature is assumed, equation d5 3b reveals that the brightness is more similar to the central 
one. Furthermore, this also applies to the projected temperature (Section ^. 3t . 

To constrain the effect of Ax, we considered the SS00 X-ray emissivity function and obtained the surface brightness profile for the set of 
objects shown in Fig. [8] This test sample well covers the temperature range in which bubbles are observed, from a to f the central temperatures 
are Tcc = {3.09, 1.83, 1.08, 0.64, 0.38, 0.18} x 10 r K. For the current analysis, it is convenient to use formula dC lit -which gives the 
spatial profile as a function of temperature- together with Table ID21 Alternatively, one can use Table ID31 to numerically integrate instead of 
summing up. For this option, the sum symbols in JcTO are replaced by \ S c #{s) with A = 3, S c #(s) = (5/2A ) # } (^ 2 - 

tf~ 1/2 )A x (tf, Z c )[(l ~ d 5/2 ) 2 ~ s 2 /i? 2 ] _1/2 dtf and A = 3 x 10~ 23 erg s" 1 . The results are shown in Fig. [9] We have that in the soft 
band the surface brightness profile evolves a follows: it begins in a limb-brightening mode for the hottest bubbles (Tcc > 10 7 K), later it 



2 It is slow except for exceptionally short times (t < 1 Myr) and low densities (no < 1 cm 3 ). 
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Figure 8. Panel (a): X-ray spatial profiles for the C95 X-ray model with Zq =Zq. The lines scale uniformly with metallicity. From top to bottom of the 
panel, the central temperatures corresponding to the lines are 3.09,1.83,1.08,0.64,0.38, and 0.18, in units of 10 7 K. If one assumes the W77 model, these 
temperatures might correspond to objects with the following parameters: lines a-d, SSCs with masses of 10 7 , 10 6 , 10 5 and 10 4 Mq an age t = 3-Myr, i.e. 
i?38 = 3000, 300, 30 and 3 (Leitherer & Heckman 1995), respectively; line e, OB association with E-^s = 0.3; and line f, W77 "typical interstellar bubble" 
blown by a star with Ess = 0.0126. For all cases no = 10-cm -3 and Vg = 1. The reverse shock position is not included in the calculation. The top arrow 
indicates the direction of the profile time evolution for a fixed set of parameters (i?38,no). The bottom arrows show how the profile varies with the energy 
input rate, the ISM density and the central temperature. Panel (b): spatial profile for cluster c in panel (a) when a constant emissivity ( a = 0, solid line) and a 
bremsstrahlung dominated one (a = 1/2, dotted line) are considered. The respective pseudo-triangle approximations (Appendix let are marked with squares. 
The reverse shock position was included in these models. 



passes through an 'unstable' almost flat mode that occurs when the central temperature is close to Tc c ~ 6.5 x 10 K (where the X-ray 
emissivity has an absolute maximum for Z >Z©), and it ends in a limb-darkening mode (centre-brightening) for lower central temperatures 
(Tbc ~ 1-2 x 10 6 K). At early ages, tc, <C 1, the limb-brightening can be extremely strong due to both an increase of the main peak nominal 
value and the visual effect produced by the relative closeness of the reverse and leading shocks. This last effect would give the central region 
a void appearance, if it were not for the contributions of regions A and B. 

Finally, we found that the surface brightness was much smaller in the hard band. For stellar blown bubbles, the hard band emission was 
practically negligible. In shape, the profiles were similar to the one shown in Fig. 1 121 



3.3 Projected temperature 

As it was indicated in Paper I, most of the superbubble X-ray emission is in the soft energy band because of the contribution of the outermost 
layers with temperatures in the range of T ~ 5 x 10 5 -10 6 K. However, the projected temperatures can be fairly larger and more similar to 
those near the centre of the superbubbles. In some cases, they can be of the order of 1 keV. Superbubbles with hot interiors (T ~ 10 7 K) 
have already been reported by Townsley et al. (2001, 2003) and Dunne et al. (2003). In Fig.[l0] we present the projected temperature curves 
corresponding to Fig. [9] It follows that according to the models, the spectra of shell-evaporating superbubbles could be more adequately 
described by a two temperature plasma model: a hot component to account for the central part and most of the bubble area, and a warm 
one to account for the outer layers. These theoretical temperature maps are independent from the speed of the geometrical evolution because 
they were derived from a ratio, see equation d39t . For the hottest bubbles (i.e. when Tc c > 10 7 K), Tx(0) is a good measure of the average 
projected temperature inasmuch as their temperature maps are quite flat almost up to the contact discontinuity. Likewise, for 'warm' bubbles 
(Tbc ~ 10 6 K), the projected central value is a good approximation up to half of the projected size. Because of this, in what follows we shall 
indistinctly refer to the average projected temperature and Tx(0) as Txc- 

For the C95 model, the temperature map can be obtained from formula d53t with a — 1 normalised by itself with a — 0. The central 
value can be obtained in the same way from equation J54b . Its low and high temperature limits are noteworthy: Txc — > (l/3)7c c as 
Tcc ^> Tcut and Txc — > T> as To — > T cut . Here, the temperature shift is an exclusive consequence of the thermodynamical structure 
established by evaporation. Our weighting scheme considers additionally the effect of Ax and predicts a central temperature (see Appendixes 
HJandO 

•9f=£cut f ^ „ \ f= min ('?i>'?Rs) 



E 



E 



> 3 + 2a , 

Tx c = — . : — . ,, - : . (55) 



v 1 + 2a , , , 

where #s and #f are the normalised to Tcc upper and lower temperature limits for which the power law with index a is valid. For a realistic 
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Figure 9. Spatial profiles in the soft band for the sample shown in Fig. [8] Here, the SS00 emissivity was considered. The contributions from hydrogen (dashed 
lines) and metals (dotted lines) were separated. The solid lines represent the sum of both contributions. The dark lines show the profiles when R B is included 
in the calculations. The red lines are obtained when the reverse shock position is neglected. Notice that by our definition, the bubbles in panels (a) and (b) are 
not well-developed (i.e. R B /R P > 0.4). None the less, they have high central temperatures; thus, the position of the reverse shock does not affect their X-ray 
luminosity significantly (see Fig. [7). 



emissivity the intervals and indexes are indicated in Table lD2l The sum starts on the interval in which Tr s falls. This means that the reverse 
shock position is embedded in the formula; none the less, given that the projected temperatures at intermediate projected radii are alike to the 
one obtained at the centre when R s is neglected (Fig.llOt. it is also useful to remove its effect. This is accomplished by setting $r s = 1 (i.e. 
R s = 0). Notice that as the coefficients X a depend on metallicity so does equation l !55t . The value of Txc for different central temperatures 
and abundances is shown in Fig.Qj] The best-fitting equations are also displayed. Note that the value of Txc decreases as the metal abundance 
increases. This occurs because the emission at lower temperatures increases with the metal abundance. 

According to our analytics (i.e. no observational artefacts) Txc just depends on the central temperature. Nevertheless, it should not be 
assumed to be To, because it is only a fraction of the later. It should not be taken either as a temperature similar to those of the outer layers. 
This is particularly important for interpreting observations of hot bubbles with Txc ~ 10 7 K (see Section [6T2t , 



4 BUBBLES WITHOUT EVAPORATION 

If the ISM encircling the superbubble owns an intrinsic strong magnetic field, the supposed dissemination of cold material to the hot bubble 
interior could be inhibited (Band & Liang 1988; Soker 1994; Dunne et al. 2003). The magnetic field prevails even after the material in 
question has been trapped in the cold outer shell and their components tangential to the contact discontinuity could critically hinder the 
evaporative process. In this case, the expression for the X-ray luminosity during the sweep-up phase is very similar to equation (10) in 
Paper I, but with a coefficient of 8.05 instead of 9.5 and the proper correction for the reverse and main shock positions with respect to the 
superbubble completely adiabatic stage. 

In this case, the integral that gives the surface brightness is elemental 

(i- ')"■-£(. ' Y1. 

V R\) Rp \ max(s,R B ) 2 J 

where 



0) = A C ,n C v(t)Ax,-23(T C ,Zc) 
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Figure 10. Projected temperatures in the soft band for the bubbles shown in Fig. [8] and Fig. [9] The hottest bubbles tend to present a more uniform profile. 
However, the ratio of their projected temperatures to the central one is lower. In spite of this, in some cases their projected temperatures can be of the order of 
1 keV. 



T XC7 =0.6270T°' 8675 , T y <3 
T XC7 =0.6452T°' 8395 , T ? ;>3 




Figure 11. Bubble projected temperature in the soft band as a function of Tq c and metallicity. Dashed line: Zq = OZ©, solid line: Zq =Lq and dotted 
line: Zq = 5Z@. The reverse shock position was neglected. Analytical expressions for the best-fitting lines (squares) to the hydrogen component and to the 
solar metallicity case are shown. Txc7 and TV are the projected and central temperatures in units of 10 7 K, respectively. As expected from the shape of the 
emissivity function, the projected temperature decreases with increasing metallicity. 



Ao^t) = 2 x 10- 23 (1.454) 2 i -^+^^-g = (1.80 x ir/Vgs" 1 pc^^K" 1 , (57) 

and we have used the Koo & McKee (1992b) average post-shock density. In this case, the average temperature of the shocked gas is close to 
the temperature at the star cluster surface, Tc ~ T sc , because a fraction of the thermal energy is used to displace the outer shell of swept-up 
interstellar material. For very fast winds, this hot plasma is expected to emit in the hard X-ray band. The corresponding surface brightness 
spatial profile, Sc.nev, is given by the quantity between square brackets in equation \56\ . It is presented in Fig. 1121 
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Figure 12. Surface brightness spatial profile for a superbubble generated by a cluster of 10 5 Mq when no evaporation occurs. Here, tg = 10 (-E38 = 30), 
no = 10-cm — 3 and Vg = 1. The solid line traces the profile considering the reverse shock position. The dashed line indicates how the inner region would 
look if the reverse shock is neglected. A slight limb brightening near the reverse shock is expected together with a smooth transition to a limb darkening 
towards the contact discontinuity. 



5 A COMPREHENSIVE X-RAY PICTURE 

In Paper I, we used the total X-ray luminosity to compare the emission of superwinds and their associated superbubbles. In this paper, the 
corresponding X-ray models have been improved by refining the analytics and incorporating explicitly a realistic emissivity function that 
separates the contributions from hydrogen and metals. Here, we give a description of how to proceed in order to obtain the theoretical X- 
ray picture. All the parameters appearing in our formulae are expected to be derived observationally. The models can handle observational 
artefacts as long as they could be parameterized as power laws or series in the temperature, as in Mazzotta et al. (2004, see Sections [2.3l and 

m. 

The theoretical diffuse X-ray luminosity of the star cluster core can be obtained from equations d32b or ( 133b . The theoretical X-ray 
luminosity of the free wind is ~ 25 per cent that of the core and its extension can be estimated from equation dA7b . The superwind X-ray 
surface brightness profile is given by equations d36b and {37} which may be helpful to extrapolate the profile of the observed diffuse emission. 
Estimators of the wind spectroscopic temperature (Tx) are given by equations d40b and d41b . All the estimators yield that for solar metallicity 
Tx is very close to the core central temperature, which is given by d21b . The last estimator depends explicitly on Z. As the polytropic index 
was left as a free parameter, these models admit a large variety of wind configurations, as explained in Section[2] 

The bubble X-ray luminosity can be obtained from d51b if the Weaver et al. (1977) model is assumed; however our approach is more 
general. For bubbles, once that the observational data would have been reduced and analysed, one of the first things to do for comparing with 
any model that considers thermal evaporation from the outer shell is to derive the theoretical central temperature, To, from the one obtained 
for the observed spectrum. This can be done using the formulae in Fig.[TT]or equation d55b . With To, the shape of the theoretical surface 
brightness profile is automatically known because it depends just on the former. For a power-law X-ray emissivity, the surface brightness 
profile can be obtained from d53b . For a more general function, equations dCl lb and dC12b can be used. Particularly, we have calculated 
the profiles that correspond to the X-ray emissivity tabulated by Strickland & Stevens (2000). These profiles are shown in Fig. [5] and can 
be used as a reference. Additionally, the derived value of To can be helpful to constrain the parameter space of the bubble model that one 
is considering, or well, to derive some parameters in order to compare them with the observations. In this regard, To, the observed X-ray 
luminosity and the model-independent X-ray luminosity given by equation d48b can be used to compare the observationally derived electron 
density with the one that follows from a shell-evaporating bubble model. All of this can help to decide whether or not a specific model is 
adequate to describe the observations. 

The X-ray luminosity of bubbles without evaporation was given in Paper I. The corresponding surface brightness is given by equations 
d56b and d57b . The temperatures of their hot interiors are similar to the temperatures at the edges of the star clusters that generate them, 
equation d 1 8b . 

The main conclusion of our previous paper was that in X-rays, the superwind contribution dominates just for compact and powerful 
SSCs while the bubble dominates in the rest of the cases. Nevertheless, that conclusion was derived from the cumulative X-ray luminosity, 
which is not directly observable, but estimated from the photon counts on the detectors. What is directly observed is the surface brightness. 
Even when the bubble could dominate by cumulative X-ray luminosity; generally, its emission would be diluted in a volume much larger 
than that of the SSC, and as a consequence, its surface brightness would be smaller than that of the superwind (as long as the objects could 
be spatially resolved). 

In Fig. Q~3] we show the combined surface brightness of a fast superwind (Vs = 1) and the bubble that it blows. Panels (a), (b) and 
(d) correspond to an evaporative bubble, and panel (c) to a bubble for which evaporation has been suppressed. In panels (a), (b) and (c), the 
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Figure 13. Combined X-ray surface brightness of regions A, B and C. In panels (a), (b) and (c) the superbubble is blown by a 1000 km s — 1 wind generated 
by a cluster with a radius of 1-pc and a mass of 10 5 Mq. The ISM density is 10-cm -3 and Z = Zq =Z,q. Panel (a) shows a shell-evaporating bubble with 
i?38 = 30 and t§ = 3 (R p » 161-pc and Tq c Ri 1.08 X 10 7 K). This bubble corresponds to line c in Fig. [9] and thus, it exhibits limb-brightening near 
R p . In panel (c) the same E^g and tg are used but evaporation is not present; hence, the bubble has the same size than in panel (a), but it is hotter and exhibits 
limb-darkening towards R p . Its average temperature is Tq 4.42 X 10 7 K. Panel (b) shows the effect of age and the smaller energy deposition rate at early 
times (Leitherer and Heckman 1995). The used parameters are E$g = 10 and t$ = 0.75. This produces a smaller bubble with R p ss 56-pc; however, the 
bubble temperature is very similar (Tq c ~ 1-06 X 10 7 K) to that of the bubble in panel (a). In panel (d), the surface brightness profile is almost flat because 



the ISM has a larger density and the star cluster is more extended and less massive: the cluster has a mass of 10 Mq and R. s , 



10, Vg = 1, te 



E38 = 3, Z = Zq =Z 



© , n ■ 



100-cm -3 , R p = 64-pc and T Cc = 7.3 X 10 6 K. 



cluster mass is 10 5 Mq, its radius is 1-pc and the ISM has a density of 10-cm -3 . Panels (a) and (b) represent the same bubble at different 
age. Its relative surface brightness (with respect to that of the wind) decreases as it expands. Since the surface brightness of the superwind 
goes as E 2 and that of an evaporative bubble goes as E^ 35 (according to the W77 model), the later effect is enhanced by the increment of 
E by a factor of ~ 3 after 3-Myr (Leitherer & Heckman 1995). Similar arguments apply for bubbles without evaporation. The parameters 
used in panel (a) are the same that were used in panel (e) of fig. 6 in Paper I, where it was shown that the bubble dominates by total X-ray 
luminosity[f| Here, we demonstrate that because of the core compactness, the superwind dominates in the X-ray picture. In panel (c), the 
bubble has the same parameters that in panel (a), but evaporation is not present. The same dynamical model is assumed; thus, the bubble 
has also the same size, but it is hotter due to the thermalisation process at the reverse shock. Finally, in panel (d) we show the interaction 
that results when the ISM has a larger density and the cluster has a larger size and a smaller mass. For the assumed parameters, the bubble 
dominates and the profile is almost flat. 

Thus, for compact (when R sc is up to several parsecs) and powerful (E^g > 1) clusters, the central emission is dominated by the cluster 
wind (region A). Around the core, there is a still bright and rapidly decreasing component due to the free wind (region B), but most of the 
brightness at radii larger than a few times R sc and up to the reverse shock comes from background emission from the superbubble (region 
C), as the contribution from the free wind becomes negligible in comparison. The rest of the projected emission (r ^ i? s ) comes from the 
superbubble itself and adopts the profiles shown in Fig.|9]as it evolves with age. 



3 The X-ray luminosity estimated in Paper I has to be corrected by a factor of ~ 2 to account for the SS00 realistic emissivity function. However, this does 
not affect the current discussion. 
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For powerful but extended, or compact but less powerful clusters, the brightness transition from one region to the next is smoother 
because the central peak produced by the wind emission becomes smaller and wider. Additionally, high ISM densities also contribute to 
smooth out the transition, since the brightness of region C depends on n 3 , 1 '' 35 (W77). Only when this circumstances are enhanced, the bubble 
surface brightness can be comparable to or even higher than that of the superwind. 

Nevertheless, there is an additional factor that one must consider. If the evaporation process is inefficient, the bubbles will show limb 
brightening, but the overall amplitude of their emissions will decrease (see Sections [3.2| and |6T2l . In such a case, the corresponding pictures 
would look similar to those shown in panels (a) and (b) of Fig.Q~3] This means that if the integration times are not large enough, the observed 
filling factors will be small, even for young normal clusters (i.e. for clusters that do not qualify as SSCs). This would translate into observing 
centre brightening or even an apparent spatial disconnection between the X-ray emission and the motion of the outer shell. 



6 COMPARISON WITH OBSERVATIONAL DATA 
6.1 Nearby clusters 

Here, we give a brief summary of the clusters analysed in Paper I. In Table[T] we give the X-ray luminosities obtained from formula ( 1331 ) for 
the same sample. The first six columns indicate the observed values. The last three give the results from our analytics. The contributions from 
hydrogen (Lxa.h) and metals (Ioca.m) were separated. The amplitude of the central surface brightness (Aab) is given in the last column. 

NGC 3603. Using color-magnitude diagrams for the massive stars contained within the core, Sung & Bessell (2004) determined an 
age of rsc ~ 1 ± 1 Myr for this cluster. Given the remarkable similarity of the spatial distribution of its stellar content with that of R136, 
the core of 30 Dor in the LMC (Moffat, Drissen & Shara 1994), we adopt a mass of 1.7 x 10 4 M@ (see Paper I). This corresponds to an 
energy deposition rate of 1.7 x 10 38 erg s _1 (Leitherer & Heckman 1995). Moffat et al. (2002) studied the X-ray emission of this cluster and 
determined a core radius of 0.71-pc. They reported a diffuse X-ray emission of 2 x 10 34 erg s _1 in the 0.5-10 keV band. Their best-fitted 
spectra have temperatures Tk c v ~ 3.1 keV for the core and kT ~ 2.1 keV for the outer region. They also reported an approximately 
gaussian-shaped X-ray surface brightness, same that is in agreement with equation {36}. The utility of having theoretical estimations for the 
cumulative X-ray luminosity and surface brightness profiles of star clusters is manifest from that work: its authors used the exterior profile 
of the diffuse X-ray emission to extrapolate to that of the cluster core, where the high density of point sources made less reliable its direct 
estimation; even after the use of removal algorithms. 

Arches Cluster. The Arches cluster is a very compact and young star cluster. Figer et al. (1999b) determined its age to be rsc ~ 2 ± 1 
Myr. Adopting a low-mass cut-off of 1-Mq, they estimated a total stellar mass of A/gc ~ 1.1 x 10 4 M Q . This corresponds to an energy 
deposition rate of Ess =1.1 (Leitherer & Heckman 1995). The cluster radius is ~ 0.2-pc (Figer et al. 1999a). Its X-ray emission was studied 
by Law & Yusef-Zadeh (2004), whom reported a total X-ray luminosity of (0.5 - 1.0) x 10 35 erg s _1 with Z = (4 - 5) Z©. Just 40 per 
cent of this is diffuse emission, being the rest due to point sources. As in Paper I, we adopt the Law & Yusef-Zadeh (2004) one-temperature 
plasma model; thus, Tk c v ~ 1.5 keV. 

Quintuplet Cluster. For the 1-pc radius Quintuplet cluster, Figer et al. (1999b) determined a mass of M sc ~ 8.8 x 10 3 M Q using again 
a low-mass cut-off of 1-M , The cluster age has been estimated to be rgc = 4 ± 1 Myr (Figer et al. 1999a). According to the Leitherer & 
Heckman (1995) evolutionary models, at this age, supernovae have already started to contribute to the wind energy budget. Hence, we adopt 
E = 2.64 x 10 38 erg s" 1 . Law & Yusef-Zadeh (2004) studied the X -ray emission of this cluster with the aid of Chandra. Their best fitted 
spectrum for the cluster core corresponds to a single-temperature solar-metallicity plasma with Tk c v = 2.4 keV. They observed an X-ray 
luminosity of ~ 1 x 10 34 erg s" 1 . 

Thanks to the additional physics synthesized in equations d32M33l>, i.e. the direct incorporation of a realistic X-ray emissivity function 
and the projected temperature, we reach a better agreement between the theoretical and observed luminosities. Their ratios are now reduced 
to factors < 2 for the case of the Arches and Quintuplet clusters. A comprehensive interpretation of the X-ray emission of these clusters 
can be found in Wang et al. (2006). For the case of NGC 3603, the cluster wind X-ray emission accounts for ~ 30 per cent of the observed 
one. This prediction agrees with the results obtained by Li et al. (2006). Their observed and simulated intensity profiles also agree with the 
general form predicted by equation i36i , although each with different indexes ai and az. The later difference is due to the exponential stellar 
distribution used in their hydrodynamical model. However, our predicted bulk X-ray luminosity is consistent with theirs within a factor of 
~ 2, since the rapid decrease of the gas density warrants that most of the X-rays come from within the star cluster core and its vicinity. Also, 
notice that, as they have pointed out, non-CIE effects become important in the supersonic wind region. 

The luminosity of the metals surpasses that of hydrogen for the case of the Arches and Quintuplet cluster, whereas for NGC 3603 they are 
comparable. The effect of compactness and a relative lower temperature can be appreciated in the case of the Arches cluster. This object has 
a radius just 3.5-5 times smaller than those of the other clusters in the sample, and its temperature is about 1—1.5 keV lower. Comparatively, 
such differences do not affect too much the X-ray luminosity, but they have an important impact on the surface brightness: the Arches Cluster 
diffuse X-ray emission is more than 2 orders of magnitude brighter in projection than those of NGC 3603 and the Quintuplet Cluster, albeit 
its total bolometric X-ray luminosity is just ~3 times larger. This later factor is mainly due to the Arches higher metallicity. 
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Table 1. Comparison with Nearby Clusters 



Object ( a > 


Tc (b) 
[keV] 


p(c) 
Use 

[pc] 


E(d) 

[10 38 ergs" 1 ] 


0(e) 

[Z ] 


-^XA.obs 

[lO^ergs- 1 ] 


j (s) 

^XA,H 

[lO^ergs- 1 ] 


j (h) 

^XA,M 

[lO^ergs- 1 ] 


A (i) 
^AB 

[10 34 erg s- 1 pc" 2 ] 


NGC 3603 


3.1 ±0.7 


0.71 


1.7 


1 


1.6 


0.22 


0.27 


0.80 


Arches 


1.5 ±0.2 


0.20 


1.1 


4 


2.0 - 4.0 


0.41 


3.90 


430 


Quintuplet 


2.4 ±0.5 


1.0 


2.6 


1 


1.0 


0.71 


1.02 


1.44 



(a): object name; (b): central temperature; (c) star cluster radius, (d): mechanical luminosity; (e) metallicity; (f): observed X-ray luminosity; (g)-(i): predicted 
X-ray luminosity contributions from hydrogen and metals and amplitude of the central surface brightness. The respective units are indicated between square 
brackets. 

6.2 Bubbles: the case of M 17, magnetically hampered thermal conduction? 

M 17 -also known as the Omega Nebula- is powered by a young open cluster with an age of ~1-Myr. The nebula has an overall diameter of 
~10— 12 pc, and coinciding with it, diffuse X-ray emission has been detected by Chandra (Townsley et al. 2003). Dunne et al. (2003) have 
studied its properties using a previous data set from a ROSAT observation. Their data sets have similar parameters. The best-fitting spectra is 
described by a single temperature plasma with T spcc ~ 0.7 keV, Z =Zq and an rms electron density of n e ~ 9 x 10 -2 cm -3 . From the 
observed bubble size (they assumed a filling factor /f = 0.5), expansion velocity, and stellar content, they derived that no «40-60 cm -3 , 
E38 = 0.1 and t§ ~ 0.13. Using these parameters, they made a comparison with the W77 model. They found that nc c ~0.2-0.4 cm -3 , 
T C c ~5.0-5.6xl0 6 Kand L XC ^(3-5)xl0 35 ergs" 1 . 

The predicted X-ray luminosity differs by a factor of 120-200 with the observed one, Lxc.obs ~ 2.5 x 10 33 erg s . To explain these 
discrepancies, Dunne et al. (2003) proposed the total inhibition of the evaporative process by the strong magnetic field (100-500 uG) that 
permeates M 17 (Brogan et al. 1999). They required the winds to be significantly clumped and homogeneously mixed in order to match the 
predicted post-shock temperature (8 x 10 7 K) with the spectroscopically derived one. An alternative explanation is that while the magnetic 
field might render the evaporation ineffective, it would not be able to completely avoid it, allowing a smaller fraction of cold gas to mix 
with the hot bubble interior. In this scenario, the extrapolated central temperature of the bubble would be more similar to the post-shock 
temperature. In fact, from our analysis, we have that for their value of T spcc the almost central temperature of the evaporative bubble is ~2 
x 10 7 K (see Fig. I lib , which makes the discrepancy with the full-evaporative case predicted Tc c even worse. 

Nevertheless, if we ignore the particular prescribed dynamical evolution and just consider a general shell-evaporating bubble, we have 
that Txc ~ Tspcc- In such a case and according to Fig. II II and equation d55t . the value of Txc derived by Dunne et al. (2003) would 
correspond to a bubble with Tcc ~ 4-5 Txc (the upper bound considers the reverse shock position). This increment of temperature with 
respect to the one derived from the W77 model has to be compensated by the same decrement in density. Since Lx oc Uq c this would reduce 
the predicted X-ray luminosity by a factor of 16-25. If additionally, we consider that the C95 X-ray model overestimates the value of Lx 
by a factor of ~2 (for our derived Tc c , see Fig. [7j with respect to the averaged SS00 emissivity, we get a correction factor ~32-50, which 
would bring Lx acceptably close to the observed value. The density near the centre would be also similar to the value derived by Dunne 
et al. (2003). This would require a decrease of the classical thermal conduction coefficient (Spitzer 1956) of about 2 orders of magnitude. 
As it was indicated in Section[4] such a drastic change is possible when strong (> luG) magnetic fields are involved (see Markevitch et al. 
2003 and references therein). Such magnetic fields can hamper the thermal conduction process and establish steeper temperature profiles. 
Note that this effect is independent from the other observationally deduced parameters (no, E38 and to) because the bubble size and internal 
pressure are not affected by the correction factors introduced above (see Section^. 

The previous analysis remarks the relevance of the projected temperature as a diagnostic tool, as well as the necessity of a broader-scope 
assessment of the thermal conductivity efficiency for the case of bubbles. 

7 SUMMARY 

Our findings and their immediate consequences are: 

(i) A new and simple star cluster wind model has been developed. Assuming a polytropic equation of state we reproduced the effect of 
radiative cooling on the hydrodynamical profiles, in consonance with previous numerical calculations. Moreover, the door has been left open 
to further applications like the parameterization of other effects, the constraining of the parameter space through threshold lines and the 
analysis of other branches of the solution. 

(ii) The superwind X-ray luminosity has been calculated using a realistic X-ray emissivity that allowed the separation of the contributions 
from hydrogen and metals. Good agreement was obtained when we compared with observational data. Also, transformation laws to include 
very general emissivity functions and completely characterize the global properties of the models have been given. 

(iii) We have calculated the corresponding surface brightness and shown that although the adiabatic superwind model predicts a very 
extended X-ray emitting volume, the X-ray bright super star cluster cores are expected to be observed surrounded by compact and faint 
haloes (if their bubbles are not considered). The weighted temperatures associated to the cluster cores were also obtained. 
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(iv) An analogue to the C95 X-ray luminosity model has been obtained considering a realistic emissivity function. The model is inde- 
pendent of the assumed dynamical evolution for bubbles with R s /R p < 0.4. Significant deviations were observed with respect to the C95 
model when To -> 10 6 K. It was also found that for under-developed bubbles (when the ratios of the reverse and leading shocks are > 0.4) 
with extrapolated central temperatures on this range, the bubble X-ray luminosity can be drastically reduced by the relative closeness of the 
reverse and main shock positions. 

(v) The X-ray surface brightness of evaporative bubbles has been obtained in terms of elementary functions by introducing a pseudo- 
triangle approximation. It has been shown that limb-brightening, flat and limb-darkening modes are all possible according to the value of the 
central temperature. Because we have separated the self-similar spatial profiles from what we call an intensity amplitude and a weighting 
emissivity factor, they are valid for all evaporative models based on the classical thermal conduction theory. 

(vi) It has been demonstrated that moderately hot, well-developed shell-evaporating bubbles (To ~ 10 K) have X-ray weighted tem- 
peratures somewhat similar to their central ones (T X c ~ 0.5-1 To) in major fractions of their projected areas. On the other hand, very hot 
bubbles (To > 10 7 K) are expected to have point- wise projected temperatures of about 0.25-0.50 To - Nevertheless, the hottest bubbles can 
still have projected temperatures of the order of 1 keV. Since Txc is a measure of the temperature at which the plasma is emitting the most, 
it is important to account for the difference between this and the actual temperature at different radii. This can be also useful to evaluate the 
efficiency of the thermal evaporation process. 

(vii) We also studied the case of when the evaporative process is totally inhibited. In this case, the bubbles have very hot interiors since 
their temperatures are similar to the post-shock temperature. Their surface brightness was also obtained. 

(viii) We have assembled the X-ray brightness models and presented a very descriptive and comprehensive X-ray picture of the interaction 
of the superwind with the ISM. 

All these tools are expected to be helpful at the time of comparing the models with the observations. 
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APPENDIX A: SUPERWINDS X-RAY LUMINOSITY 

It is convenient to initially adopt a constant X-ray emissivity Ax = Ao and express the integral J24b in terms of the velocity instead of the 
radius. The related quantities will be labeled with the subindex '0'. 

Al Initial formalism 



Substituting d24b and Q} in J 1 6b and introducing the variables x = U' /U, z = —3(2?? + 1)17 and / = (1 — zx) we get that Lxao = 
X A (h + 12) with 



Region A 




(Al) 
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h = -2(2t7 + 1)U 5/2 / A x 3/2 f~^dx, (A2) 
Jo 

and 

= Tvfu ■ (A3) 

In what follows f will be called a "base", for reasons that will be evident in brief. The above integrals can be expressed in terms of Gauss 
hypergeometric functions (see for instance, Whittaker & Watson 1996 and Lebedev 1972): 

2 F 1 (a, b; c; z) = f^^Z fe) f ~ 'T'^O- ~ ^)~ a dx, (A4) 

where T is the Gamma function. It follows that a 1 — 10/3, &i = 3/2, ci = bi + 1 and a 2 = di + 1, 62 = bi + 1, C2 = 61 + 2. Then, 



J_ Fi _ 2 £l±iW 

26i 61+I 



(A5) 



Above, to simplify the notation we have used the label of each family of indexes to make reference to the corresponding hypergeometric 
functions: Fi — 2-Fi(ai, 61; ci; z) and F2 = 2F1 (02, 62; C2; 2:). We shall keep this custom in what follows. Each U G [0, 1 / (277 + 1)] yields 
the value of Lxao at the radius indicated by l !16t . Alternatively, one may use the temperature T = T c (l — U) £ [T c , T sc ] as the independent 
variable. 

For rj — I , the total luminosity can be approximated by 



mi-) in-' 1 -1\ Ax,-23(T, Z)E% S w ln 33 -1\ Ax,-23(r, Z)E^ S 

LxAo.totai = (1-12 X 10 ergs ) ' = (9.63x10 ergs ) — '— ' , (A6) 

rtsc,pc Vg -ftsc,pc kcV 



where we evaluated dA5t at (/ = C/ S c and took Ao = A%(T, Z). The value of the emissivity function for the reference temperature T is 
taken from the SS00 tables. 



Region B 

For quasi-adiabatic superwinds, r\ « | for r > J? sc and (7 £ [t/ S c, 1). In principle, superwinds with VooP of the order of 1000 km s~ 
very extended X-ray envelopes, since the connection between temperature and radius is given by 

T \ ~ 1/4 / T x ~ v/2 



* = ¥ J { ¥c J ■ CAT) 

Nevertheless, the outermost layers meagerly contribute to the luminosity, since most of the emission from this region comes from a 
compact shell of gas around the star cluster. Thus, it suffices to estimate the external X-ray luminosity and surface brightness of quasi- 
adiabatic superwinds because those of superwinds in the radiative regime would be alike, despite the effective cropping of their X-ray 
emitting envelopes by radiative cooling. Introducing the variables x = U' /Ui and z = Ui we find that 

1 u,=u 



■3/4 



;UF a -^-F 4 



(A8) 

u,=u ac 



Lxbo — XbAoU, 

[2(64 + 1) 46 4 
with (24 = —?7/2, 64 = 3/4, C4 = 64 + 1 and 123 — 0,4 + 1, 63 = 64 + 1, C3 = 64 + 2 and 

Xb = - f 6 - (A9) 

The velocity at the cut-off radius is U cu t = (1 — T cut /T c ). For fast winds (7 C ut ~ 1. Accordingly, we can replace the hypergeometric 
functions with their asymptotic limits using the Euler beta function and the identity r(i>)r(l — v) = 7r/sin(7rv): 

, . r( c )r(c-o-6) r( c ) -n /Aim 

hm 2 F 1 (a,b;c;z) = \ ( = — ( — — -. (A10) 

T(c — o)r(c — a) F(c — a)r(a) sin(7ro) 



3 

-3/2. 



Then, for 77 = |, F3 and F4 can be substituted by 



lim^ 1 -F 4 = 3(2- 3 / 2 7r). (AUj 



Combining these equations with l |A8l l we obtain that 

£xB0 ~ -ixAO- (A12) 

This luminosity comes from a very large volume. Since radiative cooling brings the X-ray cut-off temperature closer to 7? ac , the actual 
emitting region is more compact. Consequently, the contribution of the free wind is smaller. This conclusion does not change when the SS00 
emissivity is considered, as can be shown using the transformation laws given below and consulting the numerical results of Paper I. Because 
of this, we shall focus on region A, but providing a general formalism to account for region B if it is desired. 
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A2 Transformation laws 

Before proceeding any further we shall establish laws for transforming equations dA5b and dA8b such that they can handle more general X-ray 
emissivity functions with explicit dependence on metallicity. For this reason, Ao was kept intentionally inside of the integrals dAlb and dA2b . 
These laws can be used to calculate other integral properties of the superwind such as its mass and energy within a certain volume. The laws 
also apply for bubbles. 

Let a, b' , a and A a >,Ah>, A a be arbitrary constants. The following transformations apply to dAl b and dA2b and similar integrals like 
those for region B and C: 

Law I: if Ax = A a f a then the transformation a — > a — a applies for the leading index and Ao — > A Q for each term0 
Law II: if Ax = A a f' a with /' = (1 — xz') and z' independent from z, then 2-F1 — * Fa, where Fa is the Appell hypergeometric 
function which has two bases (/, /') and two leading indexes (a, —a'): 

F A (a,-a\b-c-z,z')= J {C) C x h ~\l - xT^ 1 f~ a f a ' dx. (A13) 

r(6)r( c -6) j 

Again, Ao — > A a i applies for each term. 

Law III: if Ax = A b ix b then we need to use first the transformations b — > b + b' and c — > c+b' and then Ao — ► T(c — b)T(b)/F(c) Ay 
for each term. 

Since a, b and c, a and a' are mutually independent, we can apply the three laws simultaneously if it is required. Because all the 
hydrodynamical variables are functions of /, /' and x (or it) these laws suffice to calculate any integral property. 

In this paper, we used them to calculate the X-ray emission assuming a realistic Ax- This was done as follows: considering the more 
general case Ax = ~}2 where A a are constants that can scale with metallicity (see Table lD2l > and each a £ [—4, 4], we applied Law 

II to the quantity between square brackets in dA5b and proceeded to expand it rationally around (7 = (after dividing by T") recurring 
to admissible hypergeometric transformations (see the references in Section IaTT ). We used a (0,2) Pade approximant to make the rational 
approximation around 0. This warrants a coincidence of up to the second order with the Taylor expansion of the original function and 
provides an accurate approximation for large values of 77, since Use = 1/(277 + 1)- To account for small values of 77 and the effect of the 
exponents of the power laws, we later introduced two curvature matching filters. The next expression valid on U £ [0, 1/(277 + 1)] was 
obtained: 

r 2 i 4/10 

350(1-10*5 1 + %^^ 

^>^)= — 9(f i + U + ^) ' (A14) 

where 

Ml><x) = 70(a + 2877 + 14), (A15) 
and 

cf> 2 (ri, a) = 17a 2 + a(501 + 9527?) + 9528^(1 + 77) + 2382. (A16) 

The quantity between parenthesis on the numerator of equation dA14b is a matching filter of the form (T/T c ) b that corrects for the 
curvature introduced by a. Likewise, the quantity between square brackets is a filter of the form [1 — g(v)U b ] c that corrects for the curvature 
introduced by 77 > 0.3. The filter parameters where obtained by minimization of the relative error of the approximation. For the SS00 
emissivity function we always obtain an accuracy better than 5 per cent. In general, the same applies for all emissivity functions with 
a £ [— 4, 4] provided that 77 ^ 3/2. This ensures a good coverage of the radiative wind case. For the most extreme combination of 
parameters we are considering, i.e. for 77 w 0.3 and a ~ 4, the accuracy is ~ 15 per cent. 

From equation dA14b and Law II, it follows that the cumulative X-ray luminosity is 

Lx A (77,{a},l7) =X A (7 3/2 ^^(77,a,(7)A a r c Q , (A17) 

a 

where X\ is given by equation dA3t . 



APPENDIX B: EVAPORATIVE BUBBLES SURFACE BRIGHTNESS 

Just for convenience, we assume initially an X-ray emissivity function with a power-law dependence on the temperature, Ax = A a T a , 
where the constant A a might or might not scale with metallicity. We then express all the relevant variables of the evaporative model (Section 
in terms of the ratio of the point-wise temperature T c to the central value: 7? = T c /Tc c . From equations d46b and ( 147b one has that 

r = R p (l - tf 5/2 ), 7i c = nee??" 1 , s = R p (l - til' 2 ), s cut = R p (l - tf^ t 2 ) and R s = R(l - fi 5 ^ 2 ). The surface brightness can then 

4 In the more general case, Ax = Ax(/), the Einstein summation convention can be applied to express the function as a series in /. Similar arguments are 
valid for the other two laws. 
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be separated into a temporal amplitude, Ac', a weighted emission factor, A Q ; and a spatial profile Q(s) -between two reference points s; and 
S{- that implicitly but slowly changes with time, such that ac = Ac(t)A a @($i', #/) with 

A c (t) = 2R p n 2 Cc Ai, 



Ac* — 



and 



7 = A a ,-23-iCc, 

Al 



(Bl) 
(B2) 



r* (j 



t9 5 / 2 )(2 - 



5/2 



where Ai = 10 23 erg cm 3 s 1 . Notice that > #f because the bubble temperature gradient is negative. 
The integral above can be expressed in terms of Appell hypergeometric functions Fa 

-, t»=0f 



#2" 



tf5/2 



F A2 



F A] 



(B3) 



(B4) 



where bj = ± + §a, F A i = F A [l/2, 1/2, b i; bi + 1; tf 5/2 /tf s 5/2 , tf 5/2 /(2 - #s /2 )], and similarly, Fa 2 = F A [l/2, 1/2, bj + 1; 61 + 
2; $ 5/ ' 2 /$s / ' 2 , i? 5//2 /(2 — #s )]• For a single power law X-ray emissivity, the surface brightness spatial profile in terms of the temperature is 



crc(0, *) = A c (t)\ a O[d s , a; min(i9 s , #« a ), tf cut ]. 
In terms of the projected radius s we have that 

ac(s,t) = Ac{t)X a Sc[s,a;max.(s,R s ),Rcut]- 



(B5) 



(B6) 



In the last formula, we used the transformations indicated at the beginning of this appendix to go back to the spatial coordinate; thus, 
0(# s ...# re /) — > Sc{s ...sxref)- For the singular cases -when both arguments of the minimum or maximum function are the same- we 
define min(a, a) = max(a, a) = a. For a realistic emissivity expanded as a power series, equation JB 5 1 immediately transforms into 



A c (t) J2 A c ,e[i?s,a;min(j? i) i?R 3 ),i?f],ifs^i? s 



(B7) 



-Ac(i) A a O[i9 s , a; min(i9 s , t?f], otherwise; 



here, i?i and #f are respectively the upper and lower temperature limit for which the power law with index a is valid. 

Equations dB4t -( IB7t and their analytical properties are enough to calculate the surface brightness spatial profile. Nevertheless, we still 
need closed algebraic formulas for evaluating them swiftly and accurately. They are obtained in the next appendix. 



APPENDIX C: HYPERGEOMETRIC PSEUDO-TRIANGLES 

Consider the negative of the quantity between the big parenthesis in equation dB4t evaluated separately at the end-points while considering 
a single power-law emissivity and ignoring temporally the reverse shock position 



(CI) 



where #f is a reference temperature below # s (rf > s). For a single power law rt = Rcut- We get then a pair of functions that we express in 
terms of the variable r = flt^ 2 = (1 — s /R p ) 

*w(r;a)= ^F A1 (»,1, ^) - ^rF A2 (", 1, ^) , (C2) 

TcAr,Tf,a)= ^FAi^a.^-^^^a (C3) 

At r = rf, the two functions above take the same value, since n is the cut-off value for the emission with temperature slope a. By 
symmetry, the functions in question also take the same value at the limit r — > l - . The simplest way to prove this is to find the primitive 
function of the integrand in equation JB3b in spatial coordinates and carry out a comparison with equation dB4t . Proceeding this way, one 
can find that 

Fc,H(l";a) = F c ,c(l~,r f ;a) = 1 = — ^— . (C4) 

01 1 + /a 
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Figure CI. Hypergeometric pseudo-triangles. The hypothenuse Ten is given by equation iC2) . The function Tc c> that corresponds to the catheti, is given 
by equation {OJ. Top: T Cc = 1-68 X 10 7 K, a = (Chu et al. 1995 model). Bottom: T Cc = 6.4 X 10 6 K, a = 1/2 (Bremsstrahlung). Solid lines: actual 
hypergeometric functions. Dashed lines: pseudo-triangle approximations. 



Thus, equations JC2b and dC3t form a closed path: a pseudo-triangle. This is per se a very interesting property both from the physical 
and mathematical point of view, which complete study we reserve for a later time. Here, we give just a brief discussion that is general enough 
to solve the (super-) bubble surface brightness problem. The actual shape of the triangles depends on r f 2//5 = i?f = Ti/Tcc- For To > 10Tf 
there is practically no curvature, whereas for To ~ T the catheti are deformed (Fig. lClt . For the range of a that we are considering (see 
Table lD2t. a series expansion shows that the hypotenuse is always, for any practical purpose, a straight line. A convenient point that one 
can use to trace it is + . There, the second term in dC2b vanishes. After working with its integrand, the first term on l lC2b can be reduced to 
(l/6i)0Fr(6i + l)iF 2 [1/2, b; 6 + 1/2; r/(2- r)]/T(bi + 1/2). At the limit r -> + , the previous Gauss hypergeometric function reduces 
to unity. As a consequence J~c,n takes the value 

r(i ±2a } 



(C5) 



r(±±p + i) 

The function G(x) = F(x)/T(x + 1/2) can be approximated on the interval (—1/2, 1) by the following (1,2) Pade approximant around 



0: 



G(x) = 



T(x) 



[tt 2 +31n(4) 2 ]s + 61n(4) 



r(x+i) ^x{[tt 2 - 31n(4) 2 ]a; + 61n(4)}' 

For the case under study, when x falls out of this range, it can always be mapped back using the relation T(x + 1) 
equation of the line that traces the hypothenuse is 



Tc,n{T-a)^ 0FG (±±^) + [jJL_ _ ^ G 



■ l+2a 



(C6) 
xr(x). The 

(C7) 



The slope of this line, m, just depends on a. For a = 0, it reduces to m(0) ~ —tv/2 — — arccsc(l) = — arcsin(l _1 ). This hints that a 
secant-like approximation can be used to obtain an expression for the catheti. Because we have a closed path, a single secant approximation 
on r G [rf , 1] just recovers the expression for the hypothenuse. We consider a multi-secant approximation of the form sin(/3) ~ siii-j^jFce — 
5/(1 + 2q)]/(1 — t)} = g(r), where g(r) accounts for the variable angle j3 that the catethi subtend with the negative direction of the r-axis. 
Consistency requires the closure of the path to be conserved; thus, it is obligatory that g(jf) — sin(— m). Taking into account the catheti 
dependence on Tf/r, g(r) can be written as g(r) — h(T{/r) sin(— m), where ft(rf/r) is a function that satisfies h(l) = 1. Many of such 
functions exist, for instance, any of the form h(rf /r)/h(l). Equation dC31 > diverges as ~ Tf/r as r — *■ 0, so we use }i(ti/t) — Tf/r. This is 
equivalent to a first order approximation. It suffices for the present discussion. Thus, the expression for the catheti is 



Fc,c(t, n; a) 



jj— + (1 — r)arcsin 1^ sin(— m)J ,ifabs(m) ^ 5; 



l+2a 



+ (l-r){>- 



i(— m)l } , otherwise; 



(C8) 



where the last case separation is necessary to take the correct value of the arcsin function, specially when machine-evaluated. For a single 
power law emissivity, equations < IB4b -( IB5b and the expressions for the hypotenuse and the catheti give the next approximation for the spatial 
profile as a function of the normalised temperature 

1 &=min(0 s ,i? Rs ) 
Ui=*cat 

where the dummy variable £ is meant to be evaluated at the indicated limits as the result of an integral would and min($ s , $ft s ) brings the 
reverse shock position back into the game. In terms of the projected radius, the profile is 



e A (# s ,a;£ f ,£i) = (24 -tfsT 1 / 2 [^ +Q ^c,c(^ 5/2 ,C 5/2 ;a)]| 



(C9) 
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Sca (s, a; & , &) = 



-1/2 



ff=max(- 



(CIO) 



The expression that corresponds to a realistic emissivity (see Fig. lDll and Table lD2t parameterized as piecewise continuous power laws 
is obtained from equation dB7t 



o"c(iM) = < 



*f=*rait 

A c (t) AaeA^s.asmin^t,^),^], ifs^fl,; 

#f=*out 

-Ac(i) AuSa^s, a; min(i? s , i9f], otherwise; 

a,iSi,i> f <iJ s 



(Cll) 



or, in terms of s 



Sf— -Rcut 

Ac(t) 22 ^cxSca[s, a; max(si, R s ), Sf], ifs^_R s ; 

Sf — -Rcut 

-^c(i) A a 5cA [a, a; max(s, Si), Sf], otherwise. 

a ,Sj ,Sf >s 



(C12) 



These formulae give the surface brightness spatial profile for all evaporative bubble models. They can be also used to calculate the 
weighted temperature, as we show in Section [331 



APPENDIX D: AUXILIARY TABLES AND FIGURES 



This paper has been typeset from a TgX/ ETgX file prepared by the author. 
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Table Dl. LIST OF RELEVANT SYMBOLS 



jyiiiutji yucLiic ) 


A/T(=* nti 1 n O" 
1VJ.CCU lillg 




A tit* nmti n 

IVlCtllllllg 


indexes 'A', 'B' & 'C 


properties at regions A,B and C, respectively 


index 'X' 


a property in the X-ray band 


index 'sc' 


values at r = R sc 


index 'c' 


central values 


index 'H' 


contribution from hydrogen 


index 'M' 


contribution from metals 


index '7' 


temperature in units of 10 7 K 


n* 


II, = {Rsc, <?e, {<?m, Voap}} 


A (6) 


surface brightness amplitude 


E(3) 


mechanical luminosity 


Ess (2,5) 


E in units of 10 38 erg s _1 


F(1S) 


Gauss hypergeometric function 2F1 


F A (19) 


Appell hypergeometric function 


fc(3) 


Boltzmann constant 


KQ) 


polytropic proportionality constant 


Lx (5) 


cumulative X-ray luminosity 


M(3) 


mass deposition rate 


n(5) 


particle number density 


no (9) 


ISM particle number density 


P(3) 


pressure 


<7m (3) 


mass deposition rate per unit volume 


?e(3) 


energy deposition rate per unit volume 


r (3) 


radial coordinate 


r d (6) 


angular diameter distance 


fl(3) 


a normalised radius: R = r/R sc 


Bsc (3) 


star cluster radius 


^?sc,pc (2,5) 


star cluster radius in parsecs 


Rent (6) 


X-ray cut-off radius 


R P (9) 


principal shock position 


R B (9) 


reverse shock position 


s (6,7) 


projected radius 


5(10,21) 


surface brightness spatial profile 


t(9) 


time 


te(2,9) 


time in Myr 


T(3) 


temperature 


T x (7) 


X-ray weighted temperature 


r kc v (4,5) 


superwind central temperature in keV 


Tn s (10) 


temperature at 


T cut (6) 


X-ray cut-off temperature 


^spec (7) 


X-ray spectroscopic temperature 


u(3) 


velocity 


(7(3) 


squared velocity: u 2 


Vooa(3) 


adiabatic terminal speed 


VooP (3) 


polytropic terminal speed 


V 8 (2,5) 


VooP in units of 1000 km s _1 


V sc (3) 


star cluster volume 


X(5) 


X-ray luminosity scaling factor 


Z(5) 


metallicity 


a (5) 


X-ray emissivity temperature slope 


7(3) 


adiabatic index 


r(i8) 


Gamma function 


A (10,21) 


pseudo-triangle approximation 


6(3) 


total energy per unit mass 


«bp(3) 


region B e in the polytropic case 


»?(3) 


polytropic index 


\a (10,20) 


weighting X-ray emissivity factor 


Ax (5) 


X-ray emissivity function 


Ax,- 23 (2) 


Ax in units of 10 -23 erg s _1 cm 3 


A Q (5) 


X-ray emissivity proportionality constant (Ax = A a T a ) 


A x (9) 


evaporatively averaged Ax 


W (4) 


mean mass per particle for a ionized gas 


(5) 


mean mass per particle for a neutral gas 


P(3) 


density 


(7(6) 


X-ray surface brightness 


i? (9) 


region C normalised temperature 


#cut (9,20) 


region C normalised T cu t 


6(6) 


object subtended angle 


(20) 


surface brightness profile in terms of 1? 


r (20) 


tf5/2 


^(5,19,20) 


combination of hypergeometric functions 



The page in which each symbol appears for the first time together with pages that contain important related formulae are indicated between parenthesis. 
The list is not exhaustive, but provides the definition of all relevant symbols. A description of the most used indexes is also given. The physical properties of 
each region (as they are defined in Fig. [TJ are labeled with the subindexes A, B and C, except when no place for confusion exist. The same applies for the 
subindex X, which designates an X-ray property. 
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BAND Component Temperature Parameters Fitting Details 

(xlO 7 K) 



SOFT 


Hydrogen 


0.05-0.1 


Ao = 


0, A a = 5.12 x 10~ 27 , a = 4 


R 2 


= 


0.991, e RMS = 


= 0.0014 


(0.3-2.0 keV ) 




0.1-0.75 


Ao = 


-0.0329, A a = 0.3802 x 10~ 7 , a = 1 


R 2 


= 


0.998, e RM s = 


= 0.003049 






0.75-3.5 


A = 


0.590, A a = -18700, a = -0.69 


R 2 




0.9993, e RMS 


= 0.001659 






3.5-25 


Ao = 


0.488, A a = -4.15 x 10" 10 , a = 1 


R 2 


= 


0.9802, e RMS 


= 0.003738 




Metals 


0.05-0.1 


Ao = 


U, A a = 3.70 X 10 , a = 4 


ti 




0.995, e RM s = 


= 0.007 






0.10-0.75 


Ao = 


4.09, A a = -5630, a = -0.53 






0.987, e RM s = 


= 0.090 






0.75-25 


A = 


0.201, A a = 1.01 x 10 11 , a = -1.54 


R' 2 




0.9993, e RM s 


= 0.0176 


HARD 


















(2.0-8.0 keV ) 


Hydrogen 


0.05-0.5 


fa 
















0.5-1.1 


Ao = 


0, A a = 1.2 x 10 — 21 , a = 2.767 


R 2 




0.9961, e R MS 


= 0.0008925 






1.1-8.5 


A = 


-0.365, A Q = 2.238 x 10~ 4 , a = 0.4621 


R 2 




0.9944, e RM s 


= 0.0143 






8.5-25 


Ao = 


0, A a = 0.0318, a = 0.1667 


R 2 




0.9274, e RM s 


= 0.01013 




Metals 


0.05-0.5 


fa 
















0.5-1 


Ao = 


0, A a = 1.54 x 10 -21 , a = 2.835 


R 2 




0.9951, e RMS 


= 0.00194 






1-8.5 


Ao = 


1.208, A a = -63.1674, a = -0.2513 


R 2 




0.9951, e RMS 


= 0.00194 






8.5-25 


Ao = 


0, A a = 8.7793, a = -0.1502 


R 2 




0.9799, e RM s 


= 0.003425 


BOTH** 


Hydrogen 


>0.9 


A H = 


1.23 x 10" 4 , a H = \ 


R 2 




0.96, e R MS = 


0.04 




Metals 




Ami 

Am2 


= 4.49 X 10 17 , a M1 = -f 
= 0.79, a M2 = 


R 2 




0.98, e RM s = 


0.04 



Power laws fitted to the SS00 emissivity tables. On each temperature interval, a normalised X-ray emissivity of the form Ax,— 23 = Ao + A a T a was fitted. 
For the metals, the values of Ao and A a scale linearly with Z. The mean square errors correspond to the dimensionless Ax _23, thus they are also in units 
of 10 — 23 erg s _1 cm 3 . The curves were fitted considering piecewise continuity as a requirement. 



Table D3. Fitted X-Ray Emissivity Function - RATIONAL FUNCTIONS 



BAND 


Component 


Temperature 
(XlO 7 K) 


Fitted A x , -23 


Fitting 


Details 








SOFT 


Hydrogen 


0.05-25 


0. 2442T^ +10. 13T^ + 7. 547^-1. 0792^ +0.05112T 7 -0.00083 
T 7 i +17.91T*+21.08T3 +16. 683^+1. 063T 7 +0. 1009 


e R MS 


= 1.40 


X 


io- 


6 


(0.3-2.0 keV ) 


Metals 


0.05-25 


21.87^-1.817^+38.26^-0.7704117-0.04085 
T^+83.31T^-129.7T^+91.13T|-19.09T 7 +1.953 


e RMS 


= 3.40 


X 


io- 


6 


HARD 


Hydrogen 


0.2-0.5 


0.07841T|+0.2146T 7 -0. 1441T 7 +0.02784T 7 -0.00113T 7 -0.0001026 
T|'+0.891T 7 1 +0.2492T^+1.052T|+1.475T7 + 1.06 


e RMS 


= 3.46 


X 


io- 


6 


(2.0-8.0 keV) 




0.5-25 


0.5726T 7 +40.17T 7 -46.13T 7 +19. 07T 7 -2.812 
T*+39.6T^ +114.87^ +224.1T 7 -17.12 


e RMS 


= 1.09 


X 


io- 


4 




Metals 


0.2-0.5 


0.20847^ -0.1051T| +0.01765T 7 -0.0009808 
T^-0.16T 7 +0.3066 


e RMS 


= 3.57 


X 


io- 


6 






0.5-25 


0.3816T 7 +1.4S7T* +6.0S9T 7 -0.1862T 7 +3.877T 7 -1.478 
T 7 '-2.751T 7 l +32.31T 7 i +31.ST 7 ! -81.11T 7 +108.5 


e RMS 


= 2.30 


X 


io- 


6 



Emissivity function used for numerical calculations. For TV < 0.2 the emission on the hard band is negligible since it is several orders of magnitude 
(^ 4) smaller than the emission at the cut-off temperature in the soft band. Computationally, the evaluation of the rational approximations is faster than 
interpolating from the tables. For all cases R 2 fa 1 within the used precision. The mean square errors correspond to the dimensionless Ax _ 23 ■ to order 
to obtain this table, we associate to each temperature T; in the SS00 tables a function e; = i (Ax(Tj) — Rjk({a}, {b}; T;) 2 , where Rj^ is a (J, k) Pade 
approximant. Then we defined E = E(e\ , e, 2, ...) and solved the associated non-linear least-squares problem min \E\ \ 2 = min(J^j e 2 ) using the trust 
region method described by Byrd, Schnabel & Shultz (1988) and Dennis & Schnabel (1996). A similar process was carried out for Table lD2l In the present 
case, the values of j and k were successively increased. The most compact approximations with R 2 > 0.9 and the lowest RMS errors were selected in the 
end. 
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Figure Dl. Comparison of the SS00 X-ray emissivity tables for hydrogen (solid lines) and metals (dashed lines) with the fitted approximations (squares). Top 
panels: rational approximations. Bottom panels: power law approximations. Panels (a) and (c): emissivities in the soft band. Panels (b) and (d): emissivities in 
the hard band. 



